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Abstract 



We formulate equations for the slow time dynamics of fluid motion that 
^-j- ■ self consistently account for the effects of the variability upon the mean. The 

£C) \ time-average effects of the fluctuations introduce nonlinear dispersion that acts 

■ to spatially smooth the transport velocity of the mean flow relative to its circu- 

lation or momentum velocity, by the inversion of a Helmholtz operator whose 
. length scale corresponds to the covariance of the fluctuations. 
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1 Introduction 

We seek equations for the slow time dynamics of fluid motion that self consistently 
account for the effects of the variability upon the mean. In formulating such equa- 
tions one must choose a suitable decomposition of the flow into its rapid and slowly 
varying components and determine a strategy for applying the corresponding averag- 
ing procedure. We consider Reynolds type decompositions of either the Lagrangian 
fluid trajectory, or the Eulerian fluid velocity. These decompositions lead respec- 
tively to either the Lagrangian mean (the time average following a fluid parcel), or 
the Eulerian mean (the time average at a fixed position). Our strategy in seeking 
self consistent slow time dynamics is to apply these decompositions and their corre- 
sponding averaging procedures to Hamilton's principle for an ideal incompressible 
fluid flow. The resulting Lagrangian mean and Eulerian mean equations are obtained 
in both cases by using the same Euler-Poincare variational framework pL 0. Hence, 
these equations possess conservation laws for energy and momentum, as well as a 
Kelvin-Noether circulation theorem that establishes how the mean properties of the 
fluctuations affect the circulation of the mean flow. 

Thus, we present two formulations of the mean equations we seek: a Lagrangian 
mean theory; and an Eulerian mean theory. It turns out these theories possess a 
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certain duality. In particular, the Eulerian mean velocity appears as the momen- 
tum, or circulation velocity in the Lagrangian mean dynamics, and vice versa. The 
effect of the averaging in either case is to make the advection or transport veloc- 
ity smoother than the momentum, or circulation velocity, via the inversion of a 
Helmholtz operator that relates the two velocities. The length scale that appears 
in this Helmholtz operator is the covariance of the fluctuations, which has its own 
dynamics in each case. Thus, the temporal averaging in Hamilton's principle leads 
to a dynamical spatial filtering in the resulting equations of motion. This is the 
main point of the paper. The two formulations we present here complement each 
other and provide a flexible unified basis for further investigation and analysis of 
the effects of fluctuations on mean fluid dynamics. 

We begin by introducing new equations that describe the Lagrangian mean ef- 
fects of advected fluctuations in 3D incompressible fluid motion. The results include 
a new second moment closure model for 3D fluid turbulence. This model describes 
dynamically self- consistent interaction between a Lagrangian mean flow and a dis- 
tribution of advecting rapid (or random) fluctuations described in slow time (or 
statistically) by their Lagrangian mean covariance tensor. Its derivation combines 
the Euler-Poincare theory of fluid dynamics, the Taylor hypothesis for advection of 
the fluctuations and the Reynolds decomposition of the Lagrangian fluid trajectory. 
We also consider the effects of rotation and stratification in this model, for the sake 
of its potential geophysical applications. To help develop intuition about the so- 
lution behavior of these Lagrangian mean models we discuss ID and 2D subcases, 
as well. Then, we compare the Lagrangian mean models with their Eulerian mean 
counterparts and emphasize the duality between them. 

1.1 Lagrangian mean equations 

This Section introduces a new self-consistent dynamical model that describes the 
Lagrangian mean effects of advected fluctuations on 3D incompressible stratified 
fluid motion in a rotating frame. This model is based on the Lagrangian fluid 
description of two standard assumptions: (1) Reynolds decomposition of the La- 
grangian fluid trajectory; and (2) the Taylor hypothesis, that rapid fluctuations are 
advected by the mean flow. We substitute these two assumptions into Hamilton's 
principle for the Euler equations of an ideal incompressible fluid, and apply the 
Lagrangian mean before taking variations. Dissipation is then introduced in the 
traditional semi-empirical fashion. The results include a new second moment clo- 
sure model for 3D fluid turbulence. The latter is a development of the one point 
closure model of Chen et al. ||]- ||, based on the viscous Camassa-Holm equations 
(VCHE), also known as the Navier-Stokes alpha model. The Lagrangian mean 
Euler- Bous sines q (LMEB) model for a stratified incompressible fluid in a ro- 
tating reference frame is given by the equations, 
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— v — u x curl R(x) + Vp + gbz = v Av , with V-u = 0, (1.1) 
where — = (— + u -v), v=(l-A)u, A = V ■ • V , 
and ^ = (V- KS -V)6 with 2k s = ^<&> = <«) • Vu + Vu T •<££) . (1.2) 



These LMEB equations include the standard Euler-Boussinesq (EB) equations as 

an invariant subsystem, for which = 0. The LMEB introduce the additional 
dynamics of the covariance of the rapid fluid parcel displacement fluctuations 
into the metric of the dynamical Helmholtz operator, 1 — A. Because of the 
advective nature of the dynamics, this Helmholtz operator commutes with the 
advective time derivative, d/dt for the Lagrangian mean velocity. 

The LMEB model introduces two different fluid velocities into the average de- 
scription. The velocities in the LMEB equations ( |1.1| ) - (1.2) are u, defined to be 



the Lagrangian mean velocity; and v = (1 — A)u, found later to be the Eulerian 
mean velocity, to order o(|£| 2 ). It is also useful to think of the velocity u as the 
particle, or transport velocity, while the velocity v is the flow, or circulation velocity. 
Because of the relation v = (1 — A)u, the transport velocity u is smoother than the 
flow, or circulation velocity v, by the inversion of the dynamical Helmholtz operator, 
(1 — A). That is, in the advective time derivative d/dt = (d/dt + u -V j , we have 

u = (1 — A) _1 v. The difference u — v = Au is the Stokes mean drift velocity, 
due to the presence of the rapid fluctuations with Lagrangian mean covariance 
The interpretation of its effects will be a recurring theme in this paper. 

The symmetric tensor k$ = ^d(££)/dt is the Taylor diffusivity tensor (de- 
fined here without including any antisymmetric corrections due to rotation). The 
displacement fluctuation £ is defined by £ = — X with X = (X^) for an averaging 
process (the Lagrangian mean) taken at fixed Lagrangian label a and denoted ( • ) 
with, e.g., 

X(a,t) = (X*(a,t;w)) = lim ^ / X*(a, t ; u) duj . (1.3) 

Thus, a displacement fluctuation satisfying (£) = is given by 

£(X,t;w) = X*(a,t;w) - X(a,t) . (1.4) 
This displacement fluctuation has covariance 

<«>==<X*X*)-XX. (1.5) 

Here, X^(a, t;u) is the spatial trajectory of a fluid parcel with Lagrangian label 
a that is undergoing EB dynamics. We assume this motion depends on both a 
slow time t and a rapid (or random) time variation u. In taking the Lagrangian 
mean of this spatial trajectory, we average at constant fluid label a over its rapid 
variation in u> during a time interval T that is long (denoted limr^oc,) compared 
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to the rapid time scale for variation in u, but during which the slow time variation 
in t may be regarded as fixed. In particular, we shall assume that the result of 
the averaging operation ( • } is independent of the magnitude of the time interval T. 
Thus, in equation ( |Q| ) one may think of the Lagrangian trajectory X^(a, t;oS) as 
following the original EB dynamics in a flow regime in which a separation of time 
scales into slow (t) and fast (cj) is possible, and think of X(a, t) as following the 
approximate slow time dynamics determined by the LMEB equations. We shall see 
that p is the Lagrangian mean pressure. The other notation in ( |1.1| ) and ( |1.2| ) is 
standard for fluid dynamics: b is buoyancy; v is kinematic viscosity; g is the constant 
acceleration of gravity; and curlR = 2fi(x) is the Coriolis parameter, which may 
depend on position. The boundary conditions of the dissipative LMEB model are 

v = , u = and • h = on a fixed boundary. (1-6) 

When dissipation, rotation and stratification are absent, the LMEB model in 
(1.1) - (1.2) reduces to the ideal Lagrangian mean motion (LMM) equations: 



— v + Vp = 0, with V-u = 0, (1.7) 

where — = (— + u -v), v = (l-A)u, A = V • <&) • V , 

and ^(«) = <«)-Vu + Vu T .(«). (1.8) 



The boundary conditions for this ideal LMM model are: 

v-n = , u = , and • h = on a fixed boundary. (1. 

The LMM motion equation (|0]) may be rewritten equivalently as 




, +v-V v+ (u-v) -Vv +Vp = 0. (1.10) 
at / >> v ' 

Stokes transport 

Thus, perhaps not unexpectedly, the Stokes mean drift velocity u — v = Au con- 
tributes an additional transport term in this motion equation for the Eulerian 
mean velocity v. On the invariant manifold (££) = the LMEB and LMM equation 
sets recover their original forms. 

The ideal LMM model preserves the total kinetic energy, 

E = \J d 3 x(|u| 2 + (^W-^) =\ Jd 3 xu-v, (1.11) 

in which the Lagrangian covariance of the fluctuations couples to the gradients 
of the Lagrangian mean velocity. (Throughout the paper, we sum on repeated 
indices.) Conservation of this energy provides L 2 control on |Vu|, provided 
is bounded away from zero. In fact, for incompressible Lagrangian mean velocity 
u, the determinant det(££) is conserved on fluid parcels. Hence, the covariance 
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does remain bounded away from zero, if it is initially so. This statement still holds 
when order 0(|£| 2 ) compressibility is allowed. The ideal LMM model also conserves 
the domain integrated momentum, f vd 3 x. (With boundary conditions ( |1.9| ), this 
conserved momentum is also equal to J ud 3 x.) 

The ideal LMM model possesses a Kelvin- Noether circulation theorem 
showing how the Stokes drift velocity generates circulation of v. Namely, 



d 
dt 



dx 



7 (u) 



S(u) 



V Auj xVu J 



dS, 



(1.12) 



where the closed curve 7(11) moves with the Lagrangian mean fluid velocity u and is 
the boundary of the surface 5(u). By virtue of the operator commutation relation 



d % x d 
—A - A— 
dt dt 



d x 
dt' 



0, 



(1.13) 



we may also express the ideal LMM motion equation ( |1,7| ) in its alternative LMM 
form, 



du 

~dt 



+ u-Vu = -(l-A)- 1 Vp, V 



u 



0. 



(1.14) 



In this form, the effect of the advected fluctuations is to smooth the pressure gradient 
in an adaptive fashion depending on the velocity shear through the evolution of A. 
The divergence of equation (1.14) yields an elliptic equation for the Lagrangian 
mean pressure, p. The Kelvin-Noether circulation theorem corresponding to the 
form ( |1.14| ) of the LMM motion equation is 



d_ 

dt 



u • dx = - 

7 (u) J JS(u) L 



V x (1 - A) _1 Vp -dS, (1.15) 



which represents the circulation dynamics of u, rather than v. 

Thus, the presence of the fluctuations with Lagrangian mean covariance in 
the ideal LMM equations has five main effects, relative to the Euler equations: 

1. it smoothes the Lagrangian mean (transport) velocity u relative to the Eule- 
rian mean (momentum) velocity v; 

2. it introduces an additional transport of the Eulerian mean velocity v by the 
Stokes drift velocity u — v; 

3. it controls ||Vu||2 in the L 2 norm via energy conservation (or energy dissipa- 
tion, when viscosity is included); 

4. it creates circulation of both Eulerian mean velocity v and Lagrangian mean 
fluid velocity u around closed loops advecting with the fluid parcels; and 



5. 



it smoothes the pressure gradient in an adaptive fashion depending on the 
evolution of the fluctuation covariance driven by the velocity shear Vu. 
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Most of this paper is devoted to deriving the LMM and LMEB equations by 
approximating Taylor series expansions and averaging over the rapid time depen- 
dence in Hamilton's principle for the Euler equations of an ideal incompressible 
fluid. We average over "fast time" at fixed Lagrangian fluid label in Hamilton's 
principle for the Euler equations. This averaging over the rapid "microscopic" fluid 
motions allows us to extract an approximate Lagrangian mean Hamilton's principle 
whose Euler-Lagrange equations describe the remaining slow "macroscopic" fluid 
motions. After a discussion of these Euler-Lagrange equations from the viewpoint 
of the Lagrangian fluid parcel description, we use the Euler-Poincare theory of Holm, 
Marsden and Ratiu Q, 0, to develop and analyze the ideal LMM equations (1.7) - 



fll.q ) in the Eulerian description. We then add rotation and stratification to derive 
the ideal LMEB equations. Finally, we introduce dissipation according to the tra- 
ditional semi-empirical reasoning by requiring that the energy dissipation rate be 
negative definite. This both adds the fluctuation-dependent viscosity in the motion 
equation and introduces the dynamics for the Taylor diffusivity tensor k$ into the 
buoyancy equation. 

A second moment Lagrangian mean turbulence model is obtained by 
adding phenomenological viscosity z/Av and forcing F to the ideal LMM motion 
equation ( |1.7| ), so that, 

(— + u • V ) v + Vp = v Av + F , V-u = , (1.16) 



with viscous boundary conditions (1.6). Note that the Lagrangian mean fluctuation 



covariance appears in the dissipation operator A. In the absence of the forcing 



F, this viscous LMM turbulence model dissipates the energy E in equation (|1.11|) 
according to 



dE 
~dt 



j d 3 x tr(Vu T • • Vu) + Au • Au . (1.17) 



This strictly negative energy dissipation rate is the reason for adding viscosity with 
A, instead of using the ordinary Laplacian operator. 

The generalization of the Lagrangian mean motion model to Riemannian man- 
ifolds - to make it applicable in any coordinate system and to elucidate its intrin- 
sic geometrical structure - will be discussed elsewhere ||. We will also defer the 
derivation of the isopycnal and hydrostatic versions of this model to another time 
and place M. 



1.2 Eulerian mean equations 

We shall also present a parallel development of these Lagrangian mean results, 
but applied to Eulerian mean models. The Eulerian mean Euler- Bous sines q 
(EMEB) equations are given by 
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+ VjVU j - U x curlR(x) + VP t E t 

+gbz + I (U k ■ U ,) V(( k ( l ) E = uA E V, (1.18) 



2 



d / d 



where ^ = (J^ + tJ -Vj, V = (1-A B )U, A £ = V • < C0 B - V , 

§ = j t UO E = o, P t f t = P-i|u| 2 -^(cV) £ (u fc -uO ; 

and V- U = , where U = (U) B . (1.19) 



Here P is the Eulerian mean pressure, and v and k are constants representing 
viscosity and diffusivity, respectively. The boundary conditions for this viscous 
EMEB model are: 

V = , U = , and ( GG) E -n = on a fixed boundary. (1.20) 

These equations correspond to the same level of approximation as those for 
LMEB, but they are based on the Eulerian mean velocity U, appearing in the 
Reynolds fluid velocity decomposition, 

U(x, t; w) = U(x, t) + U'(x, t; u) . (1.21) 

The Eulerian mean averaging process at fixed position x is denoted ( • ) E with, 
e-g-, 

1 r T 

U(x,t) = (V(x,t;tj)) E = lim - / U(x, t ; u) duj . (1.22) 

T->oo T J_ T 

We shall show that the Eulerian velocity fluctuation U' is related to the Eulerian 
displacement fluctuation £ by 

= C • VU + U'(x,(;w). (1.23) 

Consequently, the Eulerian mean kinetic energy due to the velocity fluctuation sat- 
isfies 

(|Uf > £ = <CWU, fc .U,,. (1.24) 

Note that the advection of the Eulerian mean displacement fluctuation covariance 
(CC) E by the Eulerian mean velocity tJ is componentwise: d{C,C) E / dt = 0, so 
each component of this symmetric tensor is carried along with the Eulerian mean 
flow as if it were a scalar function. Moreover, we shall show that the momentum 
V appearing in the Euler-Poincare formulation of these EMEB equations is in fact 
the Lagrangian mean velocity for this theory. This duality between Lagrangian 
mean and Eulerian mean theories is another theme of the present paper. 

In the absence of dissipation, rotation and stratification, the EMEB model in 
( 1.18j ) - ( 1.19D reduces to the ideal Eulerian mean motion (EMM) equations: 
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idt + ° " V ) V + V]Vt)j + VP '°* + 2 ^ ' V (C fc C') E = , (1.25) 
where V = (1 - A £ )U , A B = V-(CC) B -V, 

V-U = 0, and + U .v)(CC) S = 0. (1.26) 



The boundary conditions for this idea/ EMM model are: 

V-n = , U = , and ( • h = on a fixed boundary. (1-27) 
The EMM motion equation (|1.25Q may be rewritten equivalently as 



'^ + V-V)V+ (V-U) X curlV + VP-^(CV) £ V(U fc -Ui) =0. 



Stokes vortex force 



(1.28) 



Thus, the Stokes mean drift velocity V — U = A E \J contributes an additional 
vortex force in this motion equation for the Lagrangian mean fluid velocity V. 

On the invariant manifold V( = the EMM equation set reduces to the 

3D ideal Camassa-Holm model ]l]], The Eulerian mean system ( |L25|) - (|1,26 ) 



recovers the 3D ideal Camassa-Holm model as an invariant subsystem for (C k C l ) E = 
a 2 5 kl , with a a constant length scale. 

The ideal EMM model preserves the total kinetic energy, 

E =\j d 3 x(\U\ 2 + (( k ( l ) E U >k ■ U, ( ) = l - j d 3 x U • V , (1.29) 

in which the Eulerian covariance of the fluctuations couples to the gradients of the 
Eulerian mean velocity. Conservation of this energy provides 1? control on |VU|, 
provided (££) E is bounded away from zero, which it will be, if it is initially so. The 
ideal EMM model also conserves the domain integrated momentum, f V d^x. (With 
boundary conditions (1.27), this conserved momentum is also equal to f Ud 3 x.) 



The ideal EMM model possesses a Kelvin- Noether circulation theorem 

showing how the Stokes drift velocity generates circulation of V. Namely, 



d ' V-dx 1 



dt 7 7 (u) 2 j js(v) 



V(U, fc -U,,) XV(CV} E -dS, (1.30) 



where the closed curve 7(U) moves with the Eulerian mean fluid velocity U and is 
the boundary of the surface •S'(U). 

Thus, the presence of the fluctuations with Eulerian mean covariance ( CC) E m 
the ideal EMM equations has four main effects, relative to the Euler equations: 

1. it smoothes the Eulerian mean (transport) velocity U relative to the La- 
grangian mean (momentum) velocity V; 



2. it introduces an additional vortex force into the Lagrangian mean velocity 
equation due to the Stokes drift velocity U — V = A E Xj; 
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3. it controls || VU||2 in the L 2 norm via energy conservation (or energy dissipa- 
tion, when viscosity is included); and 

4. it creates circulation of Lagrangian mean velocity V around closed loops ad- 
verting with the Eulerian mean velocity U. 

A second moment Eulerian mean turbulence model is obtained by adding 
phenomenological viscosity uA E Y and forcing F to the ideal EMM motion equation 
flLgp , so that, 



(l + U .V)V + ^W + VPf ot + ^(U, fc .U,,)V(C fe C , ) S (1-31) 

= z^A^V + F, where V-U = 0, 



with viscous boundary conditions ( 1.20| ). Note that the Eulerian mean fluctuation 
covariance appears in the dissipation operator A E . In the absence of the forcing 



F, this viscous EMM turbulence model dissipates the energy E in equation (1.29) 
according to 



dE 

— = - v I d*x 
dt 



j d 3 x tr(VU T • ( CC) E -VU) + A E U • A E U . (1.32) 



This energy dissipation law justifies adding viscosity with A E , instead of using the 
ordinary Laplacian operator. 

1.3 Outline of the paper 



In Section [2], we shall recall the Euler-Poincare equations (2.11) in the context of the 



ideal Euler fluid equations. In Section |3[ we shall introduce into Hamilton's principle 
for Euler's equations the Reynolds decomposition of a Lagrangian fluid trajectory as 
the sum of its mean and fluctuating parts. In the remainder of Section ||, we shall: 
substitute this decomposition into the Lagrangian for Euler's equations; transform 
the resulting rapid (or random) Lagrangian L(uj) in equation fl3.4| ) into the Eulerian 
description; make an approximation of it based on truncating a Taylor expansion; 
take its Lagrangian mean, denoted (L) in equation ( |3,34| ); and finally use the Euler- 
Poincare theorem to obtain approximate Lagrangian mean motion equations at order 
0(|£| 2 ). In Section || we shall provide the physical interpretations of the quantities 
arising during this modeling procedure and compare the resulting equations with 
other models. In Section |H| we shall discuss both ideal and viscous one point closure 
models based on the Euler-Poincare equation. We shall also compare these one 
point closure models with Reynolds-averaged Navier-Stokes (RANS) models, Large 
Eddy Simulation (LES) models, and the Leray regularization of the Navier-Stokes 
equations. 

In Section |6] we shall augment the Euler-Poincare equation in ( |6.71 ) to include 
dynamically varying correlations of the rapid/random fluctuations and use this La- 
grangian mean fluid theory to derive a new approximate motion equation for the 
slow evolution of the Lagrangian mean fluid velocity in Eulerian coordinates. This 
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is the basis for the Eulerian analysis of the Lagrangian mean motion (LMM) model. 
The LMM equations in (|6.17| ) - (|6.2C| ) will include self-consistent forces caused by 
the correlations of the rapid/random fluctuations and expressed in terms of their 
Lagrangian mean covariance tensor, The approximate dynamics of the co- 

variance tensor itself will also be determined, thereby producing a two point, 
or second moment, Lagrangian mean closure of Euler's ideal incompressible fluid 
equations. Several properties of this ideal second moment closure model are also de- 
rived in Section |(| including its conservation laws for momentum and energy, and its 
Kelvin-Noether circulation theorem, which are all inherited from the Euler-Poincare 
formulation. We note that the Lagrangian mean effects of the rapid fluctuations are 
purely dispersive at this stage and, thus, are energy conserving. 

The Lie-Poisson Hamiltonian structure of the ideal LMM model is given in Sec- 
tion 0. Its 2D behavior is discussed briefly in Section ||. Upon adding viscosity 
in Section ^, we shall introduce a second moment Lagrangian mean closure of the 
Navier-Stokes equations corresponding to the LMM model and compare this closure 
with the Reynolds averaged Navier-Stokes equations. 

In Section [ji], we shall discuss geophysical applications of these ideas and use 
our approach via the Euler-Poincare theory to derive the Lagrangian mean Euler- 
Boussinesq (LMEB) equation set ( |10.2| ) - (10.3) for the Lagrangian mean motion of 
a rotating stratified turbulent incompressible fluid. Section [ll] considers lower di- 
mensional subcases in one and two dimensions, in an effort to help develop intuition 
about the solution behavior of these Lagrangian mean models. 

Section 12 studies the Eulerian mean counterpart of the LMM model and empha- 
sizes the duality and parallel mathematical structures shared in the two approaches 
through their formulations as Euler-Poincare equations. Section [l3| adds rotation 
and stratification to the Eulerian mean model, then discusses some of its aspects in 
fewer dimensions. Section 14 contains a summary and conclusions. 



2 Review of Hamilton's principle and Euler-Poincare 
equations for ideal fluids 

Consider the Lagrangian L comprised of fluid kinetic energy with volume preserva- 
tion imposed by a Lagrange multiplier P (the pressure), 

+ P(X(a, t), t) ( det(X^)(a, t) - l) j . (2.1) 

Here, X(a, t) is the Lagrangian fluid trajectory: that is, x = X(a, t) is the current 
position of the material point starting at initial position a at time t = 0. We denote 
the derivatives of the function X(a, t) by X = dX/dt and X^ = <9X/9a. After a 
brief calculation, Hamilton's principle, 



L(X,X) 



d 3 a 



X(M) 



6 dtL = 0, (2.2) 
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for the Lagrangian L(X, X) yields the following Euler- Lagrange equations, 



5L 5L , . 

,5X oX. 



given explicitly by (dropping boundary and endpoint contributions), 

BP 

<5X: X + det^)— =0, (2.4) 
SP : det(X„) = 1 . (2.5) 

These are Euler's equations for the incompressible motion of an ideal fluid in the 
Lagrangian description. 

We obtain the Eulerian description of this motion by defining the Eulerian fluid 
velocity U(x, t) and volume element D(x,t) via the basic kinematic relations, 

U(x, t) = X(a, t) and D(x,t) = ( detX^(a,i)) 1 at x = X(a,t). (2.6) 



The volume element D satisfies the continuity equation, 

^ + V-( J DU) = 0. (2.7) 

The volume element D(x,t) is an advected quantity, (sometimes called "frozen- 
in"). By this, we mean a quantity that is expressible purely in terms of Lagrangian 
labels; so that it is invariant along the Lagrangian mean fluid trajectory, or equiva- 
lently, satisfies a certain Lie-derivative relation For example, the volume element 
Z?(x, t) is expressible as the Jacobian for the transformation from Lagrangian to Eu- 
lerian coordinates. Thus, the volume element satisfies 

Dd 3 x = d 3 a, (2.8) 

and the continuity equation ( |2.7| ) is implied by the corresponding invariance re- 
lation for such an advected quantity, in this case, 



U/V) ~ ( ¥t + £u ) {Dd " x) = (^bW + v ' (jDU) ) ur V) ' ' 2S 



where £\j denotes Lie derivative with respect to the vector field U(x, t), the Eu- 
lerian fluid velocity. The basic relations ( |2.6[) allow one to transform the Lagrangian 
( |2.1| ) into Eulerian variables as 

L = j d 3 x|^|U(x,t)| 2 + P(x,t)(l- J D(x,t))| . (2.10) 

Hamilton's principle 5jdtL = for a Lagrangian L(U, D) defined in this way 
yields the following Euler- Poincare equation for U, see Holm, Marsden and 
Ratiu 0, [|, 
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d \ 1 SL 1 SL , SL , 



(throughout, we sum on repeated indices) and the Lagrange multiplier P imposes 
the constraint, 

§-»■ < 2 - 12 > 



Substituting the variational derivatives of the Lagrangian ( 2.10 ) 



1 SL 6L 1, l9 SL , 



into the Euler-Poincare equation fl2.1l|) with constraint ( 2.12| ) and using continuity 



(2.7) yields Euler's equations for the incompressible motion of an ideal fluid in the 

+ U-V) U + VP = 0, V-U = 0. (2.14) 



Eulerian description, namely 

d_ 

dt 

For Lagrangians of the type ( |2.1[ ) — specifically those Lagrangians that are in- 
variant under "particle relabeling" (by the right-action of the volume preserving 
diffeomorphism group acting on the tangent space of the fluid parcel trajectories) 
- the Euler-Poincare equations ( |2.11| ) are equivalent to the Euler-Lagrange equa- 
tions (|2.3|) . See Holm, Marsden and Ratiu [|J], Q for more details, discussions and 
proofs of this type of equivalence. 

Of course, the Euler equations ( 2.14[) could also be found by directly trans- 



forming the Euler-Lagrange equations ( |2.4| ) from the Lagrangian, to the Eulerian 
description by using the basic relations (|2.6| ). This would avoid the step of trans- 
forming Hamilton's principle into the Eulerian description. And this is the point: 
the equivalence of the Euler-Lagrange equations (^^) and the Euler-Poincare equa- 
tions ( p. 11 ) facilitates easy, immediate transitions from one description to the other 



that are helpful in developing approximate fluid models and interpreting their solu- 
tion properties. In this paper, we shall take both the Lagrangian, and the Eulerian 
viewpoints, switching from one to the other whenever it facilitates our purpose in 
the development of these models. In the end, we shall be writing Eulerian represen- 
tations of approximate fluid models, found by taking either Lagrangian, or Eulerian 
means of Hamilton's principle and then applying the Euler-Poincare theory. 

3 Averaged Lagrangians and Euler-Lagrange equations 
3.1 Lagrangian fluid trajectory fluctuations 

The trajectory of a Lagrangian fluid parcel X*(a, t ; u>) may be decomposed into its 
mean and fluctuating parts as 

X«(a, t;u) = X(a, t) + £(X(a, t),t ; u) . (3.1) 
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This is the Reynolds decomposition of a Lagrangian fluid trajectory. Here £ = 
— X is a rapid (or perhaps random) vector field of fluctuations defined along the 
Lagrangian fluid trajectory. The independent variable ui in equation ( |3.1[ ) denotes 
either rapid time variation, or random fluctuations. So the variable u is allowed a 
dual interpretation. In its first interpretation, u is regarded as a short time scale 
associated with rapid fluctuations. In its second interpretation, co may be regarded 
as a random variable associated with a stochastic process obeyed by the fluctuations 
and defined along the Lagrangian fluid trajectory X with label a and slow time 
variation t. This variability may be regarded equally well as being either intrinsic, 
or extrinsic. We may assume the the fluid trajectory is fluctuating intrinsically; 
or, we may assume the fluid is subjected to a random forcing, which introduces a 
stochastic component into its acceleration, resulting eventually in the decomposition 



(3.1) of the Lagrangian fluid trajectory. 

We denote by (•) the averaging procedure performed at constant Lagrangian 
label a and time t. This averaging is performed either over the rapid time uj, 



(/) = lim - 



f(iv) dw , 



(3.2) 



or, alternatively, with respect to a certain probability distribution V(uj) associated 
with the random event uj, 



f(u))V{u) doj , with / V(u) doj = 1 



(3.3) 



Such an average taken while holding Lagrangian labels a fixed is called a La- 
grangian mean. The fluctuations are assumed to have zero Lagrangian mean, 
(£) = 0; so t ne quantity X(a, t) in equation is the mean Lagrangian fluid tra- 
jectory, since (X^(a, t ; uj)) = X(a, t). Thus, the quantity £ describes the fluctuating 
displacement of the fluid trajectory, relative to its Lagrangian mean. In the pres- 
ence of these rapid (or random) fluctuations, the Lagrangian L in equation (|2.1| 



appearing in Hamilton's principle (2.2) contains both slowly varying mean variables 
and rapidly fluctuating, or random variables. Thus, the Lagrangian (|2.1|) becomes 



L(uj) 



d 6 a 



1 



+ P(X*(a,t;a;),t) det(X|')(M ;uj) - 1 , (3.4) 



where superscript £ denotes the Reynolds decomposed Lagrangian fluid trajectory 
in equation (pT 



3.2 Induced Eulerian velocity fluctuations 

The decomposition X^ in (|3.1|) of the Lagrangian fluid trajectory into its mean and 
fluctuating components implies a corresponding decomposition of the associated 
Eulerian velocity field. Via the basic kinematic relations for the fluid velocity, 

U(x + £(x,«;cj),i) = X € (a,t;w) for x + £(x,t ; uj) = X*(a,t ; u) , (3.5) 
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we have 

U(x + £(x,t;w),t) = ^(x + £(x,t;a;)) = u (x, t) + ^(x, t ; w), (3.6) 

where — = — + u(x,t)-V, (3.7) 

defines the advective time derivative. Since the displacement fluctuations of the 
fluid trajectory have zero Lagrangian mean, (£) = 0, the corresponding velocity 
fluctuations will also have zero Lagrangian mean, so that (d£/dt) = 0, as well. 
Therefore, upon taking the Lagrangian mean of equation fl3.6|) we find the relation 

u(x,t) = (U(x + £,*)) = X(a,t). (3.8) 



Thus, the quantity u (x, t) appearing equations in ( |3.6j ) and ( |3.7|) is the Lagrangian 
mean fluid velocity at the position x = X(a, t) along the fluid trajectory with 
Lagrangian label a. 

We use the basic relations in (|3.5| ) to transform variables in the Lagrangian ( |3.4D 
from labels a, to position along the mean fluid trajectory x = X(a, t), yielding 

L(u) = J Dd 3 x!^\U(x + £(x,t;oj),t)\ 2 

+ P(x + £ (x, t;u),t)[ det(x£) det(X^) - l] J , (3.9) 

where the volume element D satisfies det(X^) = D _1 (x, t). Consequently, the time 
dependence of the product of determinants splits into two factors, 

det(x|- / ) det(X^) = D' 1 (x, t) det (J + V£(x, t ; w)) . (3.10) 
slow fast (or random) 

Hence we may rewrite equation ([0|) for the Lagrangian L(u) in the Eulerian de- 
scription equivalently as 




3.3 A Taylor series approximation to order 0(|£| 2 ) 

Until this point in the development, we have made no approximations. Now, though, 
we shall assume that the magnitude |£| of the rapid fluctuations is small enough to 
allow meaningful Taylor expansions of the unapproximated Eulerian velocity U and 
pressure P as, 

U(x + £(x,t;cu),t) = u(x,t) + ^ (3.12) 
= U(x,t) +£(x,i;w) -VU(x,t) +0(|£| 2 ). (3.13) 



Lagrangian and Eulerian mean fluctuation effects 



D. D. Holm 



17 



Taking averages and setting (£) = and (d£/dt) = in the two expressions for 
U(x + £, t) in (|3.12|) and (jjTl3|) gives 



(U(x + £(x, t ; w), t)) = u(x, t) = U(x, t) + 0(||| 



(3.14) 



where u(x, t) is again the Lagrangian mean velocity. Thus, the slowly varying 
velocities u(x, t) and U(x, t) differ only at order 0(|£| 2 ). 
A similar Taylor expansion result holds for pressure, 



P(x + £(x, t ; u),t) = P(x, t) + £(x, i ; w) -V P(x, i) + 0(|£| 



(3.15) 



Taking averages then gives a relation for the Lagrangian mean pressure, p, 

(P(x + £(x, t ; w), 0) = P(x, t) = P(x, i) + 0(|£| 2 ) . (3.16) 

Hence, the pressures p(x, i) and P(x, i) also differ at order 0(|£| 2 ). We shall discuss 
the physical meaning of equations ( |3,14| ) and ( 3.16 ) further and model their order 
0(|£| 2 ) terms in detail later, in Section |3]. Here, we only emphasize that the present 
approach ascribes all of the rapid variation in velocity U(x + t) and pressure 
P(x+^, t) to the Lagrangian trajectory fluctuation displacement £(x, t ; u); whereas 
there could be an additional Eulerian fluctuation in the functions U and P, them- 
selves, e.g., one could take U(x + £, t ; u). We shall return to the issue of properly 
representing the velocity and pressure fluctuations in Section ||, when we compare 
their Lagrangian and Eulerian representations. 

Taking the difference between equations ( |3.12 ) and ( 3.13; ) and using equation 
( 3.14| ) now yields an equation for the modulational dynamics of £(x, t;u)) in slow 
time, valid to order 0(|£| 2 ), 



§ = § + u(x , t) .V« = { -Vu(x, 4 ) 



(3.17) 



This equation is also equivalent to the vector field commutation relation, 





d_ 

Idt 



lot 



+ u(x,t)-V , £-V 



0. 



Hence, we find the remarkable relation 

d 



ih £-VA(x,t)=£-V-A(x,t), 



(3.18) 



(3.19) 



for any vector A(x, t), provided relation ( |3.17| ) holds. We shall see in a moment 
that this is an exact relation, provided the fluctuation field £ in the decomposition 
( |3.1| ) satisfies £(X(a, t), t ; u) = (£(a;cj) • ^)X(a, t), in which the dependences on 
fast and slow time variables are separated (factored). 
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3.4 Remarks about advected quantities and Taylor's hypothesis 



We shall show that, geometrically, equation ( 3.17 ) means that the slow evolution, or 
modulation, of the fluctuation vector field £(x, t; to) is invariant under the flow of the 
Lagrangian mean velocity, u(x, t). In other words, this vector field is an advected 
quantity, "frozen" into the flow of the Lagrangian mean velocity. This situation is 
analogous to either the vorticity stretching equation for the 3D incompressible ideal 
Euler equations, or the ratio of magnetic field intensity to mass density in ideal 
magnetohydrodynamics, Q. Following the solution attributed to Cauchy || fo r the 
3D vorticity stretching equation, the solution of the advection condition ( 3,17j ) may 
be expressed in components in separated form as 

ax** (a, t) 



e(X(8L,t;u),t)=Q{^ r(a;o;), with F\ = pQ 
slow fast 



da A 



(3.20) 



To verify this solution, one may substitute it into d£ A (a;uj)/dt\ a = 0, where £ y 
{ F ~ l )tC and {F' 1 )^ is the matrix inverse of F\, so that (F _1 )fF^ = 5 A and 

preservation of fluid volume 



B 

8\. The existence of these matrix inverse relations is guaranteed by 



det(X^) = det(fl) = detF = 1 



We also note the relation, 



dw 



1 v ' "' dxi 
Rearranging equation ( [3 .201) implies 

d 



r(x,t), for x = X(a, t) 







d_ 

da 



£ • "5 — f° r 



x 



X(a, t). 



(3.21) 



(3.22) 



(3.23) 



where £(a;u;) may be taken as the initial value in slow time of the fluctuation 
£(x.,t;oj). Thus, along the Lagrangian mean fluid trajectory the vector field £ • V 
with parameter oj may be expressed solely in terms of the Lagrangian labels. This 
means £ • V is an advected quantity, invariant in slow time along the Lagrangian 
mean fluid trajectory, just as indicated by the commutation relation ( |3.18| ). Fur- 
thermore, to order 0(|£| 2 ) the decomposition fl3.1| ) may be expressed using relation 

flOap as 

d 



X 5 (a,t 



X + ^(X,t;w) = l + e A (a;a; 



dc 



X(a,t) 



(3.24) 



Conversely, we may show that advection of the Lagrangian vector field £(a ; ^)'-§g 
implies the advection condition ( |3.17| ) for the slow time evolution of the Eulerian 
vector field £ • V, since (in analogy to the calculation ( |2.9D for the Lagrangian 
invar iance of the volume element), 







0_ 

dt 



£(a;w) 



d_ 
da 



d_ 

ot 
dt 



u-V^-^-Vu 



(3.25) 
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where £ u denotes Lie derivative with respect to the Lagrangian mean fluid velocity, 
u(x, t). Vanishing of the term in brackets recovers the Eulerian invariance equa- 
tion ( |mD , which we now understand is equivalent to assuming the separated form 
£(X(a, t),t ;u) = (£(a;u/) • ^-)X(a, t) in the decomposition ( |3.lD , 



Remark. The advection of the fluctuating vector field £ • V expressed in relation 
( 3.17| ), or its equivalent commutator form ( 3.25| ) is an extended form of the Taylor 
hypothesis, i.e., that small rapid fluctuation fields are swept downstream with the 
mean flow ||, flpfl . Relation ( 3,17| ) extends the traditional interpretation of Taylor's 
hypothesis to apply to advection of the rapid fluctuations of a vector quantity by 
the Lagrangian mean flow. The Taylor hypothesis is usually interpreted |jTo| , |p"l| 
as meaning that turbulent fluctuations satisfy d/dt = —Ud/dx for a sufficiently 
short time interval, where U is the (constant) Eulerian mean flow velocity in the re- 
direction. The extended form of Taylor's hypothesis we introduce here asserts that 
the advection relation ( |3.17 ) holds along each mean Lagrangian fluid trajectory and 
accounts for the vector nature of the displacement fluctuation, According to its 
derivation in the steps leading to equation (|3.17 ), this extended Taylor's hypothesis 
is valid to order 0(|£| 2 ). When u is replaced by (f7, 0, 0) with constant U the 



advection relation (3.17) reduces to the traditional form of Taylor's hypothesis, 
which is usually applied to scalar quantities, 111 



3.5 Averaged approximate Lagrangians 



After the approximations leading to the Taylor hypothesis relation ( 3.17| ) with its 
geometrical interpretation ( 3.25| ), our Lagrangian L(uj) in (3.11) becomes 



L(u) = Jd 3 x Id i|(l + £(x,f;a,).V)u(x,t)| 2 

[l + €(x > t;w)-V)p(x ) t)l[det(/ + VO-£>(x,t) \ , (3.26) 



+ 



in which p is the Lagrangian mean pressure. Using relation ( 3.23|) transforms this 
formula to the Lagrangian picture as, cf. equations ( |2.1| ) and ( |3.4| ), 

2 



L{u) 



d A a 



l + r(a;w] 



d 

da A 



X(a,t) 



(3.27) 



+ 



d 
da A 



p(X(a,t),t) 



det(XiO(a,t;w)-l 



Therefore, we may take the Lagrangian mean of this approximate form of L{oj) to 
find the following averaged approximate Lagrangian (using ) =0) 



(L) 



d A al - 



1 r.v. 



(3.28) 



+ p(X(a,t),t) 



iU ' { (I+ ^)) det(X * )_1 
d 



+ det(X^((edet(l+^)^^ ),>(X(a./)./) 
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In the Eulerian description, this averaged approximate Lagrangian becomes 

D 



(L) 



d 3 x 



|u| 2 + (C¥)(u fc -u,i 



(3.29) 



+ p 



<det( J + V£))-D + ({£ det(J + V£)> • V p }, 



where partial spatial derivatives are denoted by subscript comma, e.g., u & = du/dx k 
<9fcU, for k = 1,2,3. In preparation for using the Euler-Poincare theorem, we com- 
pute the following the following variational derivatives of the averaged approximate 
Lagrangian, 



1 S{L) 
D Su 

5D 

8p 



u-±(d k D(et l )di)u, 

.k • u .l I = - r h>l- 

(det(J + V£)) - D - V • (£ det(J + V£)> 



-p+I|u| 2 + i(^}(«A •"/) = -P/. 



(3.30) 



The natural boundary condition u • n = at a fixed boundary with unit normal 
vector n is imposed on the Lagrangian mean velocity in the course of deriving 
the Euler-Poincare equations. In addition, when taking the variations in ( 3.30| ), 
we assume in integrating by parts that the fluctuations do not penetrate a fixed 
boundary, so that £-n = must be satisfied at such a boundary. Consequently, 
we have fik{£, k £ l ) = at fixed boundaries, as well. We shall discuss variations with 
respect to after making further approximations in (L). 



3.6 Further approximations 

Stationarity of the averaged approximate Lagrangian (L) in (fr29b and the last 



equation in ( 3.30 ) imply 

D = (det(J + V£)> - V • (£ det(J + V£)> . (3.31) 
We shall approximate this relation using the identity, 

(det(J + V£)) = 1 + V-<£> + \ V- (£V • i - i • V£> + <det(V£)> . (3.32) 

Since we assume |£| to be small (at least typically, or in an average sense), the rela- 
tion ( 3.31| ) can be given a simple and concise approximate form by neglecting terms 



of order 0(|£| ), or higher, involving (det(V£)) and V-(£ det(V£)}. Namely, 



D = 1 - \ V- [<£ • VO + <|(V • 0)] = 1 - ^<^>,fci • (3-33) 



This order 0([|| 2 ) compressibility due to the presence of the fluctuations is consis- 
tent with the findings of generalized Lagrangian mean (GLM) theory, ]12| , [13], [14] 
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Upon substituting the order 0(|£| 2 ) approximations yielding relation (|3.33|) into 



the Lagrangian ( 3.28 ), we find the following averaged approximate Lagrangian, 



(L) = J d*x< 


f [iu i 2 + (u k . u ,)] + P [i - d - \{ee), k ] ' 


>. (3.34) 


This approximate 
and (3.28), 


(L) returns to the Lagrangian picture as, cf. equations 


0,(11 



(L) 



d A a 



|x ( a,t)| + <eV>(^-^) 



+ p(X(a,t),t) det(X^)-l 



(3.35) 



-det(X^p) 



<9 2 p 



da A da B da A da B 



d 2 * 1 fF -isC dp 



da c 



The two equivalent forms ( |3.34 ) and ( p. 35 ) of the averaged approximate La- 
grangian possess order 0(|£j 2 ) corrections involving higher order derivatives of the 
velocity, or pressure, that are contracted using the metric provided by the La- 
grangian mean covariance of the fluctuations, either in its Eulerian form, (£, k £, 1 ), 
or in its Lagrangian form, (£, A !; B ). According to relation ( |3.20| ), these two represen- 
tations of the covariance are related by 



<^f B )(a) = (^)(X,t)(F- 1 )^(F- 1 )f . 

This relation and the Piola identity, 





da A 



det F(F 



-i\ A 



0. 



(3.36) 



(3.37) 



provide the transformation laws required in converting between the two expressions 
( 3.34| ) and ( |3.35| ) for the averaged approximate Lagrangian. 



Remark on the covariance determinant. Equation ( p.36| ) for the transforma- 
tion of the covariance tensor as it advects along a Lagrangian fluid trajectory 
also implies the following relation for its determinant, 



det <«>(a) 



det (&)(X,t) 
(det F) 2 



(3.38) 



Consequently, the product D 2 det is preserved along fluid trajectories, 



|(D 2 det <«»=(). 



(3.39) 



Therefore, the covariance of the fluctuations cannot vanish in the course of the 
motion, provided it is initially nonzero. Thus, one may regard the principle axes of 
the symmetric tensor as describing an ellipsoid, that is carried along with each 
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fluid parcel and represents an additional "internal" degree of freedom for the fluid, 
associated with the Lagrangian mean covariance of the advecting fluctuations. This 
covariance ellipsoid may deform and change orientation as it follows the course of 
the fluid motion, but its volume (the determinant det (££)) cannot either vanish, or 
become unbounded, as long as the fluid density D is finite. 



Eulerian representation of the covariance dynamics. The time derivative 



of the relation (3.36) and the basic relation (3.22) imply 



-I 



<e^)(a) = f ^-(ee) - <^'K- - ^(^kf-^-t . (3.40) 



dt 



Hence, the Eulerian dynamics for the Lagrangian mean covariance (£, k £ l ) may be 
expressed as 




or, in vector form, 



^<«> = <«>" Vu + Vu T .<£0 . (3.42) 
at 



3.7 The Euler-Lagrange equations with order 0(|£| 2 ) compressibil 
ity 

The Euler-Lagrange equations of the averaged approximate Lagrangian ( |3.35| ) , 





dt 



5(L) S(L) 

-^--^ = 0, 3.43 



are given by (dropping endpoint contributions), 

We shall defer our discussion of <5(L)/<5X until later. Stationarity of (L) in (|3.35|) 
under variations in p implies order 0([£[ 2 ) compressibility, 

det F = 1 + \ ( det F^p)) AB + \ ^ ^Mf ^ (<P£ B >*>) • (3-45) 
Equivalently, by using equations ( |3.36| ) and ( |3.37| ) we find the expected relation, 



i=D+i(ee), k i- (3.46) 

To write the Eulerian form of the motion equation ( |3.44j ) we must transform the 
operator, 
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from the Lagrangian, to the Eulerian description. For this, we use the metric trans- 
formation formula ( 3.36 ), the definition ( |3.21 ) and the Piola identity ( 3.37| ) to find 
the transformation law for A, 



A = 



d 
da A 



A?B\ 



d 
da7 



detF 



' k da A detF 



d 
da B 



D~ 



dx k ^ ^ ^ dx l D 



(3.48) 



dx k det F dx l 

Consequently, we also have the remarkable commutation relation between the oper- 
ators d/dt and Ajj, 



-i 



d_ 

dx h 



D 



-i 



0_ 

dx 



d 









dt' 


I dt 



0. (3.49) 



This relation allows us to transform the acceleration term in equation (3.44) from 
the Lagrangian, to the Eulerian description, as 



+ u -V 1-Ad u 



or, equivalently, 



d_ 

dt 



1 A • » -V )U 



5X 



5JJ± 
SX 



o, 



0. 



(3.50) 



(3.51) 



In both forms of this motion equation, the divergence of the Lagrangian mean 
velocity is given by 



Vu 



1 dD 
D~dT 



1 



d 



Ml 



ld_ 
2~dl 



(£ h t l ),ki + o(|l| 



(3.52) 



2 - m l ),kl dt 

In its first form ( |3.50| ) the presence of the operator (1 — Ajj) in this motion equation 
smoothes the Lagrangian mean transport velocity u relative to the momentum, or 
circulation velocity v = (1 — A/))u. In its second form ( |3.51| ) the operator (1 — Ajj) 
acts to smooth the generalized force 5(L)/5X in an adaptive fashion depending on 
the Lagrangian mean covariance which in turn depends on the velocity shear 

Vu. At the moment, the generalized force 5{L)/5X is still a Lagrangian quantity. 
We shall see in Section ^ how to complete the expression of these equations by 
computing 5{L)/5X in purely Eulerian form in terms of a fluctuation stress tensor. 



Since the averaged Lagrangian in equation ( |3.35| ) has no explicit time depen- 
dence, Noether's theorem implies conservation of energy, 

dX dX 1 

da A da B J ' 

2 



E 



- / d 6 a 



X(a,t) +(i A i B )- 



d 3 a 



X(a, t )- + < ? * { ')(X, t )(F-')^(F-')f^.^}, 



d 3 x D 
d 3 x ( I 



M 2 + u fe -iu 
1 



(3.53) 
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where we have used the advection solution ( 3.36| ) for the Lagrangian mean covariance 
in transforming between the first and second lines, and then transformed into the 
Eulerian description in the last line, cf. equation ( 3.34| ). The conserved energy 
E in (|3.53| ) is the total kinetic energy of the Lagrangian mean motion equation 



( 3.44 ), including the mean contribution from the covariance of the rapid (or random) 
fluctuations. We note that E contains a term of order 0(|£| 4 ). 



3.8 Restoring incompressibility in the order 0(|£| 2 ) model 

We now simplify matters by restoring incompressibility of the Lagrangian mean fluid 
velocity, thereby producing a simplified averaged approximate Lagrangian, cf. 
equation ( |3.34 ) 



U 



+ p 



1-D 



(3.54) 



Here, to regain some simplicity in the analysis below, we dropped terms in the 
pressure constraint (U|) of order 0(|£| 2 ), but kept terms of the same order in the 
kinetic energy piece of the mean Lagrangian in ( p.29 ). The approximate (L) in 
( 3.54 ) so obtained returns to the Lagrangian picture as, cf. equations fl2.1| ), ( |3.4| ), 
jplSq ) and ( g3p , 



d A a 



/ax ax \ 



+ p(X(a,t),t)[det(X^ 



(3.55) 



The particular averaged approximate Lagrangian (L) in ( |3.35| ) also could have been 



obtained by substituting the approximate decomposition in equation (3.24) directly 



into the kinetic energy term of the original Lagrangian (2.1) and then averaging. 



Relative to the starting Lagrangian (2.1) for Euler's incompressible fluid equa- 
tions, our procedure of decomposition, approximation and averaging, followed by 
restoring incompressibility of u has merely introduced an additional term into the 
kinetic energy of the averaged approximate Lagrangian (L) in ( 3.35| ). This additional 
term involves higher order derivatives of the velocity, that are contracted using the 
metric defined by the pullback of the advected Lagrangian mean covariance of the 
fluctuations, namely, {^ A ^ B ), as given in equation ( 3.36 ). 



3.9 The Euler-Lagrange equations for the LMM model 

The Euler-Lagrange equations (^^) for the simplified averaged approximate La- 
grangian (|3.55| ) are immediately obtained as 
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and stationarity under variations mp implies volume preservation, det(X^) = 1. To 
write the Eulerian form of this equation we must transform the operator A as done 



in equation (3.48). When volume preservation is imposed, this gives 



A 



d 
da A 



AfB\ 



re 



d_ 



dx k dx l 



A 



D 



for D = 1 . 



(3.57) 



This relation allows us to transform equation ( 3.44 ) from the Lagrangian, to the 
Eulerian description, as 



+ u -V ) ( I - A ) u -!- V p - () . \x\\ h V- u - U . 



or, equivalently, 







1-A)( — + u-V)u + Vp=0, with V-u = 0. 



(3.58) 



(3.59) 



These are two equivalent forms of the Lagrangian mean motion (LMM) equa- 
tion. In its first form, the presence of the Helmholtz operator (1 — A) in this motion 
equation smoothes the Lagrangian mean transport velocity u relative to the momen- 
tum, or circulation velocity v = (1 — A)u. In its second form, the Helmholtz operator 
(1 — A) acts to smooth the pressure gradient in an adaptive fashion depending on 
the covariance which in turn depends on the velocity shear, Vu. 

Since the averaged Lagrangian in equation ( |3.55| ) has no explicit time depen- 
dence, Noether's theorem implies conservation of energy, cf. equation ( 3.53| ) 



E 



1 



1 



d 6 a 



2 ld " 



1 



Xfa,£ 



+ («)(X, ( ,(F-^(^)fg.||}. 



- / d x D 



and D = 1 for the incompressible LMM case. 



(3.60) 



Physical interpretations of u and v as the Lagrangian 
and Eulerian mean fluid velocities 



In this section we use the extended Taylor hypothesis ( |3.17 ) to interpret the velocities 



u and v in the Lagrangian mean motion (LMM) model by comparing the results of 
Section ||| with the results of the traditional Reynolds fluid velocity decomposition. 

The interpretations are as follows: u = (U) L is defined as the Lagrangian mean 
fluid velocity (denoted as (U) L in this section only) and v = (1 — Ad)u is approxi- 
mately the Eulerian mean fluid velocity (U) E . Here, A^j defined in equation (3.48) 



is the Laplacian operator whose metric is given by the Lagrangian mean displace- 
ment covariance for the advected fluctuations. The difference u — v = (U) s 
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is called the Stokes mean drift velocity. We shall show that (U) satisfies the 
closure relation (U) s = A^u. Thus, we shall identify v = u — A^u as the Eulerian 
mean fluid velocity (XJ) E = (U) L - (U) 5 . 

This section is meant to be self contained. The notation we use here for averaged 
quantities such as (U) L , (V) E and (U) s is intended to be as self-explanatory as 
possible. 



4.1 Stokes mean drift closure relation 

The Reynolds decomposition of the Lagrangian fluid trajectory introduced in 
equation ( |3.1| ) implies relation (|3.6| ) for the Lagrangian mean fluid velocity u, de- 
noted in this section as (U) L , 

U(x + £(x,t;w),t) = (U) L (x,t) + U / (x,t;w), with U £ = ^(x,t;w). (4.1) 
Note that the slow-time advective derivative, 

| = | + <U)^*).V, (4.2) 

acts only on (x, t) dependence and does not act on u dependence. As before, the 
angle brackets as in (U) L denote average over dynamical behavior in u that is either 
rapid in time, or random. In comparison, the traditional Reynolds fluid velocity 
decomposition is expressed at a given position x in terms of the Eulerian mean 
fluid velocity, (U) E as 

U(x, t; u) = (U) E (x, t) + U'(x, t; u) . (4.3) 

We note that the rapid w-dynamics appears in different functional forms on the 
left hand sides of equations (4.1) and ([4 .3j) in the definitions of (U) L and (U) E . 



The Reynolds fluid velocity decomposition in ( [4 .3; ) describes fluctuations in velocity 
without reference to a Lagrangian parcel trajectory. In principle, such a trajec- 
tory could be obtained by integrating this decomposition from an initial reference 
configuration. Rather than developing the analysis relating these velocity formulas 
from a Lagrangian viewpoint, we shall instead follow [^] in presenting an heuris- 
tic argument relating these velocities as Eulerian quantities by using Taylor series 
approximations and asymptotics in the magnitude |£| of the fluctuations. 
Taylor expansion of the Reynolds fluid velocity decomposition fl4.3j ) yields 



U(x + £t;w) = U(x, i ; + | • V U(x, t ; W ) + ^U ifei (x, t ; W ) + 0(|^| 3 ) 
= <U) s (x,t)+£. V(U) s (x,t)+ i^(U)^(x,t) 

3\ 



U'(x, t ; w) + £ ■ V U'(x, t;u) + 0(|£| 3 ) . (4.4) 



Upon assuming we may equate the different functional forms of velocity arising in 
the two decompositions as U(x + £(x, t; u>),t) = U(x + t ; u), comparison of the 
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two formulas (|4.1|) and (4.4) gives 



0(1) 
0(\£\ 2 ) 



(U) L 

dt 
(U) s 



(U) E ; (4.5) 

£ • V<U) £ + U'(x,i; W ); (4.6) 
(U) L - (V) E 

«-VU / >(x,0+ ^ k Z l )(U) E ki(x,t). (4.7) 



There is some risk in making this order by order comparison of Taylor expansions, 
because of the possibility of double counting the fluctuations when comparing their 
Eulerian and Lagrangian representations. However, there are two cases of the order 



0(|£[) relation (4.6) that are unequivocal: 



1. All the fluctuation is modeled in the Lagrangian fluid trajectory. (This is the 
case studied here.) Then U' = and (|4.6| ) implies d^/dt = £ • V u. This is 



the Taylor hypothesis in (3.17). 



2. All the fluctuation is modeled in the Eulerian fluid velocity. Then U = 
d^/dt = and we obtain U' = — £ • V(U)' E from equation ( [4. 6|) . The rela- 
tion d£/dt = is the Taylor-like hypothesis introduced in []23[| . We shall 



investigate this situation in Section 12, when we deal with Eulerian mean fluid 
theories. 

For the Stokes mean drift velocity (U) s = (U) L — (U) E in the present situation 
with purely Lagrangian fluctuations (so that U' = 0) we shall follow |15| and take 



(U) 5 (x,i) 



'X- vuO 

'v-UO -v)(u) L 



i 



- «v-00-v(u) L + ^(ez l )(u) L k i + o(\t\ 3 ) 

= ( V- < eO -V ) (U) L + o(|£| 2 ) , (Stokes drift) . 
We drop the latter two terms in this expansion, by arguing that 

<V.£)=o(|£|) and (U) L fc/ = o(l). 



(4.8) 



(4.9) 



Hence, we arrive at the approximate closure relation for (U) s , to order o(|£| 2 ) 




where the operator A is defined by A = (V-(££)-V) = Ad + 0(|£| 4 ), so we need 
not distinguish between A and Ad here. 
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Identification of v as the Eulerian mean velocity. Subject to the assumptions 
(|4.9| ) and accepting the definition of (U) s as in (4.8), the momentum, or circulation 
velocity v may be identified as the Eulerian mean velocity (\J) E , by following the 
chain of relations, 

(XJ) E = (U) L - (U) 5 = (1 - A D )(U) L = (1 - A D )u = v. (4.11) 

This establishes the physical interpretation of the quantity v as the Eulerian mean 
velocity. 



Lagrangian mean, and Eulerian mean pressures. A similar Taylor expansion 
may be performed to determine the relation between the Lagrangian mean, and the 
Eulerian mean pressures. Namely, 

P(x + £) = P(x) + £ • VP(x) + \ (t k H l )P >kl (x.) + 0(|£| 3 ) . (4.12) 

Hence, we obtain 

<P) s (x) = <P) L (x) - I(^)(P)^ (x ) + 0(|£| 3 ), (4.13) 

which relates the Eulerian mean pressure (P) E to the La grangian mean pressure 
{P) L to order 0(|£| 3 ). 



Remarks on boundary conditions. With the interpretation given in ( PH ), the 
Eulerian mean velocity v = (U) E naturally must be tangent to fixed boundaries, so 
that 



v • n = , at a fixed boundary. (4-14) 

The Lagrangian mean velocity u = (U) L must satisfy boundary conditions that 
allow determination of u from v by inverting the dynamical Helmholtz operator, 
1 — A. For this, we shall choose Dirichlet boundary conditions on the Lagrangian 
mean velocity, 

u = , at the boundary. (4-15) 

We note that n • ( = at the boundary, as well, since the fluctuations should 
not penetrate the boundary. These remarks provide the rational for choosing the 
boundary conditions ( |1.S| ) for the ideal LMM model. 



Relation to the Monin-Yaglom formula. An alternative formula for the Stokes 
mean drift velocity due to Monin and Yaglom [17] introduces the diffusivity tensor 
k, according to 



(U) s = V • K 



(4.16) 
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Consistency of this formula with relation (|4.1C| ) allows one to identify the diffusivity 
tensor ft as 



* = <a>-v<u) J 



(4.17) 



We decompose this tensor into its symmetric and antisymmetric components, as 
k = ks + ka, where, 



K S 
K A 



<«> -v (u) L + (v(u; 



L\T 



<«> -V (U) L -(V(U) L ) T .(^)). 



We note from equation ( [3.42 ) that 



and by direct manipulation we find 



(4.18) 
(4.19) 

(4.20) 



V . KA = --curl <£x {£■ V)(U) L ). 



(4.21) 



Relation ( 4.20h for k$ is precisely the definition of the Taylor diffusivity jl8| , ]l9[ 
and relation ( 4.21 ) implies that ka does not contribute to the divergence of (U) . 
Thus, we find 



v. ( u)« = v.. s = v.(v.i|(«>). 



(4.22) 



In combination with the LMM equations, these identifications provide dynamical 
equations for the diffusivity in the Monin-Yaglom formula which are of potential use 
in modeling geophysical fluid dynamics |2( 



Comparison with the Craik-Leibovich equations [15], [16]. The explicit 
expression ( 4.10 ) for the Stokes mean drift velocity (U) s = u — v allows us to write 



the LMM motion equation ( 3.58|) in the notation of this Section as 



-|<U) B + (U) E -V (XJ) E + (U) 5 -V (U) E + V(P) L 



with (XJ) S = A(U) L , and V- (U) i 



0. 



0, (4.23) 
(4.24) 



Thus, the Stokes mean drift velocity introduces an additional transport of Eule- 



rian mean velocity into the ideal LMM motion equation ( 4.23 ). The corresponding 
Kelvin-Noether circulation theorem is 



d_ 

It 



(XJ) E ■ dx 



V{U)f xV (U)§ -dS, 



(4.25) 



for any surface ^((U)^) whose boundary is the closed curve 7((U) L ) moving with the 
Lagrangian mean fluid velocity (U) L . Thus, circulation of Eulerian mean velocity 
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will be generated according to these equations whenever the gradients of the Stokes 
drift velocity and Eulerian mean velocity are not colinear. In this circumstance, the 
Stokes drift velocity generates convective circulation of the Eulerian mean velocity 
that tends to reduce the gradient of the Lagrangian mean velocity. 

The motion equation ( 4.23j ) of the LMM model is reminiscent of (but different 



from) the ideal Craik-Leibovich (CL) equations [15], [16]. In the CL theory, the 
rapidly oscillating waves at a free surface are assumed to be unaffected by the more 
slowly changing currents below. The effect of the waves on the Eulerian mean 
velocity is parameterized in the CL theory by introducing into the Euler equations 
a "vortex force," expressed in terms of a prescribed Stokes drift velocity (U) (x, t). 
The CL equations are given by, 

^(U} £ + (U) E -V (V) E - (XJ) S x curl (XJ) E + Vw = , (4.26) 
V-(U) B = 0, and m = P+^\(XJ) E + (U) s \ 2 -^\(U) E \ 2 . 

Here w is a modified pressure term that includes the Eulerian mean pressure P as 
well as the increase of the kinetic energy of the fluid due to the waves. The term 
(U) x curl (U) E is the "vortex force" of the CL theory of Langmuir circulation. The 
Eulerian mean fluid velocity (U) is assumed to be divergenceless and is required 
to be tangential to fixed boundaries of the domain of flow. 

The first difference from the CL equations is that the ideal LMM model has 
a dynamical equation for the Stokes mean drift velocity, while this is a prescribed 
function in the Craik-Leibovich model. Also, the Stokes mean drift velocity appears 
as a transport term ((tr) 5 -V (U) £ ) in the LMM motion equation (4.23), while it 



appears as a "vortex force" (— (U) s x curl (U) E ) in the CL model. Both of these are 
substantial differences between the LMM model and the Craik-Leibovich theory. 

4.2 An order 0(|£| 2 ) model with div{XJ) E = 

Rather than the divergenceless condition V- (U) L = appearing in the LMM equa- 
tions (p..7|), the condition V- (U) = is implied by the Eulerian mean of the 
original divergence free condition V-U = 0, when combined with the velocity de- 
composition ( [4 .3| ) and the assumption that the Eulerian mean velocity fluctuation 
has zero Eulerian mean, (\J') E = 0, where {-) E denotes average over the rapid 
dependence at fixed Eulerian position. 

If this were so, then preservation of the condition V- (U) E = would determine 
the pressure p in the LMM motion equation ( 4,23| ) by solving the Poisson equation 



— Ap = V- ((U) -V (U) s ) with Neumann boundary conditions obtained by taking 
the normal component of the motion equation ( 4.23j ) and using (U) E ■ n = at 



the boundary. This situation is appealing because of the clear physical interpreta- 
tion of the velocity variables and their boundary conditions. However, it requires 
V- (U) L = 0(|£| 2 ). Hence, to explore this possibility we must begin by restor- 
ing the 0(|£| 2 ) compressibility induced by the fluctuations and calculated earlier in 
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equation (|3.33|) as 



D = 1 



(4.27) 



The arrangements in Hamilton's principle needed for the pressure to impose this 
relation as a dynamically consistent constraint are discussed in Section || leading to 
the Lagrangian in equation ( |3.34 ) . It would be convenient if these arrangements were 
also to imply that V- (U) E = o(|£| 2 ) ~ 0, i.e., that the Eulerian mean velocity were 
incompressible to a certain approximation. We shall now investigate this possibility. 

In the notation of the present section, the continuity equation for D as in equa- 
tion (fL27l) implies, to order 0(|£| 4 ) 

Combining this relation with equation (|4. 10 ) and the divergence of equation (|4.11|) 
gives a formula for the divergence of the Eulerian mean velocity, 

V-(U) S = V-(U) L - V-(U) 5 = V-(U) L -V-(V-(£0-V(U) L ) 



2dt 



1 

2 V - 



dt 



(ur-v, v- v- 



0(|VV(U> L |)0(|^ 



(4.29) 



Thus, V-(U) £ ' = o(|£| 2 ) ~ 0, provided |VV(U) L | = o(l), as we argued earlier in 
equation ( |4.9| ). Since the Lagrangian mean velocity (U) L is smoothed by Helmholtz 
inversion, this assumption may be plausible in certain flow regimes. If this asump- 
tion is made, the approximate equations that result are 



where 



and 



^-(U)'--V : 1 +V/> = (). V-(U)^(). 



(U) i = (l-A)- 1 (U) £ , A = (v.<«>-v) 



^ + (u) L -v ) <«> = ( -V (U) L + ( V(U) L ) T • ( «> . 



The boundary conditions for the system ( |4.30| ) - ( [4.31 ) are 
(U) £ - n = , (U) L = , and • n = on a fixed boundary. 



(4.30) 
(4.31) 
(4.32) 
(4.33) 



The energy for this system is, 

E= l -fd*x{l- \{et)jk) (U) £ • (U) L . (4.34) 

This energy will be conserved exactly for V-(U) £ ' satisfying equation Q4.2EQ without 
approximation, as we found in the second order model in Section |||. And it will be 
conserved approximately for V-(U} £ ' = 0. 
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Summary of Lagrangian mean models with order 0( |£[ 2 ) compressibil- 
ity. There are two contending theories with order 0(|£| 2 ) compressibility: one has 
exactly divergenceless Eulerian mean velocity (U) and only approximate energy 
conservation; and the other has small divergence of (U) E and exact energy conser- 
vation. In fact, the conserved energy in the latter case contains a term of order 
0(|£| 4 ), which of course is higher order than the validity of the equations. The 
choice between these two theories probably should be made on the basis of their 
performance in practice. The divergenceless model defined in equations ( 4.30| ) - 
(|4.33j ) may have the advantage of being easier to implement numerically than the 
nondivergent theory, which consists of equations ( 4.30| ) - (4.33) with V-(U) £ ' ~ 
in equation ( 4.30|) replaced by equation ( 4.29 ). 

Under mild conditions on the viscous version of the model in equations 

(|4.30| ) - (|4.33j ) will dissipate the energy ( 4.34 ), if the motion equation is modified to 
introduce viscosity, as 



d_ 

dt 



+ (U) L -V )(XJ) E + Vp = vA(XJ) E , V-(U) 



E 



0. 



(4.35) 



If needed, an additional dissipative modification of the Lagrangian mean covariance 
dynamics could also be proposed. Namely, we propose 



(! + (u} L .v)<a> = UO-v<u> L + (v<u> L ) T -<tt> 



T 



(4.36) 



In this equation, r is a relaxation time, a is a length scale below which the effects 
of fluctuations on the mean flow should be suppressed and A is a diffusivity that 
suppresses gradients of ( The quantities t, a and A may all be taken as constant 
parameters, with perhaps A ~ a 2 /r. For r _1 > |V(U) L |, the additional dissipation 
terms in equation Q4.36| ) will cause the covariance ( to approach the isotropic, 
homogeneous conditions represented by the VCHE model, which will be discussed 
briefly in Section |5| At boundaries one may take to satisfy equation (4.36) 

with the A term absent. 



5 One point closure equations 

5.1 Euler-Poincare equation for the approximate (L) 



The Euler-Poincare equation ( [2.11 ) can be written for the approximate (L) in equa- 
tion fl3.54|) for incompressible flow as 



0_ 

dt 



v + (u.V)v + D j VM J + VP foi = 0, V-u = 0, 



where the momentum conjugate to the velocity u is given by 



1 6{L) 



D 5u 



D=l 



u- (0* )u. 



(5.1) 



(5.2) 
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The total pressure Ptot in equation (5.1) is defined as, cf. equation (|3.30|), 



^^-^H'-^^Xufu,,). (5.3) 

In this Section, the fluctuation covariance is taken as being independently 

prescribed, and thus is not varied in Hamilton's principle. Equations ( |5.1| ) with 
definitions Q5.2Q and ( |5.3D are generalizations of the n-dimensional Camassa-Holm 
(CH) equations derived in [||, Q. The latter equations are recovered when the 
isotropy conditions 

(Z k £ l ) = a 2 8 kl , (5.4) 
hold and, moreover, the statistics are homogeneous, so that a 2 is constant. As we 



shall see in Section 12, there is a natural extension of the CH model for Eulerian 
mean fluid dynamics, in which u is the Eulerian mean fluid velocity and v is its 
Lagrangian mean counterpart. 



Background. The derivation jy], [g], of the CH equation (5.1) with definitions 



(5.2) and (|5.3[) , and homogeneous isotropic Eulerian fluctuation statistics satisfying 



(|5.4| ) with constant a 2 generalizes a one dimensional integrable nonlinear dispersive 
shallow water model [ |21j] , [^] , to the n-dimensional situation and provides the inter- 
pretation of a as the typical Eulerian mean amplitude of the fluctuations. See [23| 



for the extension of that derivation to Riemannian manifolds and discussions of al- 



ternative boundary conditions for the case of homogeneous statistics. See also |24| 
for further analysis of the Riemannian case of the n-dimensional CH equations. 

Holm, Marsden and Ratiu |jj , [§] note that the conditions of isotropy (^J) and 
homogeneity (constant a 2 ) need to be modified near fixed boundaries, due to the 
physical requirement that = be satisfied, where n is the unit vector normal to 
the boundary. However, the condition 

™fc(£ fc ^) =0 on the boundary, (5.5) 
implied by this physical requirement, cannot be satisfied for constant a 2 in equation 



(5.4). Chen et al. |3|]- |5[, overcame this difficulty by allowing spatial variation of 



a near fixed boundaries in straight pipe and channel geometries. 

5.2 Relation to one point turbulence closure models 

We note that the velocity v defined in equation (^^) is the momentum conjugate 
to the velocity u (or dual to u, in the sense of variational derivative of the kinetic 
energy). On that basis, Chen et al. ||- [||, proposed the following viscous variant 
of (|5.l|), in which viscosity acts to diffuse the momentum v 



d 

— v + (u • V)v + Vj Vw> = i/Av - VPtot , V-u = 0. (5.6) 
at 
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In their case, Chen et al. defined the momentum v and the modified pressure P by 

v = u- (V-(O)u- {d k a 2 8 k )vi, (5.7) 

and 

Ptot =P~ ^|u| 2 - -^{ u ,k -u fc ) . (5.8) 

Chen et al. (3|- allowed for spatial variation of a 2 , particularly in flow regions 
near boundaries. They also allowed for the mean (£) to be nonzero near boundaries 
to account for the anisotropy there. Chen et al. referred to equation ( |5.6| ) with defi- 
nition ( |5.7| ) as the viscous Camassa-Holm equations ( VCHE) , although this model is 
also known as the Navier-Stokes alpha model because it reduces to the Navier- 
Stokes equations when alpha is absent. They proposed a one point closure model for 



turbulent flows in pipes and channels by comparing equation (5.6) with the Reynolds 
averaged Navier-Stokes equations in those geometries and identifying corresponding 
terms. They then verified the predictions of this closure model by comparison with 
experimental data at high Reynolds numbers. 

Chen et al. ||- Q, also gave a continuum mechanical interpretation to their 
VCHE or NS-a closure model, by rewriting (|5.6|) (in the case where the isotropy 



conditions Q5.4j ) hold, with a 2 = constant) in the equivalent constitutive form, 

^ = divT, T= -pl+ 2v(l - a 2 A)D + 2a 2 D , (5.9) 

with V u = 0, D = (1/2) (Vu + Vu T ), ft = (1/2) (Vu - Vu T ), and co-rotational 
(Jaumann) derivative given by D = dD/dt + DCl — CID, with d/dt = d/dt + u-V. 
In this form, one recognizes the constitutive form of VCHE or NS-a as a variant of 
the rate-dependent incompressible homogeneous fluid of second grade p5| , p6| , 
whose viscous dissipation, however, is modified by the Helmholtz operator (1 — 
a 2 A). There is a tradition at least since Rivlin [27] of modeling turbulence by using 
continuum mechanics principles such as objectivity and material frame indifference 
(see also For example, this sort of approach is taken in deriving Reynolds stress 



algebraic equation models l29fl . Rate-dependent closure models of mean turbulence 
such as the VCHE or NS-a closure model have also been obtained by the two-scale 
DIA approach []30| and by the renormalization group methods [31]. We shall see 
in Section ^| that the covariance contributes a fluctuation shear stress term in 
the total stress tensor for the order 0(|£| 2 ) compressible Lagrangian mean motion 
(LMM) model. 

5.3 Comparison of VCHE or NS-a with LES and RANS models. 

Reynolds-averaged Navier-Stokes (RANS) models of turbulence are part of the clas- 
sic theoretical development of the subject |l(J, p2| , [32jj. The related Large Eddy 
Simulation (LES) turbulence modeling approach |34] 



, 1 35], [gg, provides an op- 
erational definition of the intuitive idea of Eulerian resolved scales of motion in 
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turbulent flow. In this approach a filtering function ^"(r) is introduced and the 
Eulerian velocity field is filtered in an integral sense, as 

u(r)= f A'f(r-r')U B (r'). (5.10) 

This convolution of Ug with T defines the large scale, resolved, or filtered velocity, 
u. The corresponding small scale, or subgrid scale velocity, u', is then defined as 
the difference, 

u'(r)=U B (r)-u(r). (5.11) 

When this filtering operation is applied to the Navier-Stokes system, the following 
dynamical equation is obtained for the filtered velocity, u, cf. equation (|5i 



d — 

— u + u-Vu = -divT- Vp + i/Au, V-u = 0, (5.12) 

in which p is the filtered pressure field (required to maintain V • u = 0) and the 
tensor difference 



T=(V E V E )-uu, (5.13) 

represents the subgrid scale stress due to the turbulent eddies. This subgrid scale 
stress tensor appears in the same form as the Reynolds stress tensor obtained from 
Reynolds averaging the Navier-Stokes equation. 

The results of Chen et al. §§— ||, may be given either an LES, or RANS inter- 
pretation simply by comparing the constitutive form of the VCHE or NS-a closure 
model in ( |5.9Q term by term with equation ( |5.12 ), provided one may ignore the dif- 



ference between Eulerian mean, and Lagrangian mean velocities as being of higher 
order. Additional LES interpretations, discussions and numerical results for forced- 
turbulence simulations of the VCHE model will be presented elsewhere |57|. 

5.4 Comparison of VCHE or NS-a to Leray's equation 

The Leray regularization of the Navier-Stokes equations is given by p^] . 

^U + (U) r VU = i/AU - Vp, V-U = = V-(U)/, (5.14) 
in which the velocity field U is transported by the spatially filtered velocity 

(U)/(r)=r 3 / sfrHr-rOJUOrOdV. (5.15) 
Jk 3 



Here g € C°°(Br) is a smooth positive function that vanishes outside a finite sphere 
and is normalized to unity, J R3 g = 1. Thus, the spatial scales in (U)e smaller than 
t have been smoothly filtered out. The spatially filtered velocity (U)^ satisfies the 
important inequality 

|<U)<(r)| <max|U(rO|, (5.16) 
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obtained by taking the sup norm of its definition ( |5.15| ) and using the normalization 
of g. Solutions of the Leray equation fl5.14Q satisfy various regularity properties 
that are also shared by solutions of the VCHE []39| . The only difference between the 
forms of these two equations is that the VCHE or NS-a closure model contains the 
additional term Vj'Vu 3 , stemming from its derivation as an Euler-Poincare equation. 
This additional term ensures the Kelvin- Noether circulation theorem for the 
VCHE. That is, the VCHE or NS-a closure model satisfies 

d f , f r<9v • 

— h U • V V + Vj\U J 

dt 



dt 



v • dx 

7(u) J-liv) 



dx. = v<h Avdx, (5.17) 
/ 7 (u) 



for any closed curve 7(11) that moves with the Eulerian mean fluid velocity u. In 
comparison, the Leray equation satisfies 



d_ 

dt 



<p U-dx=(f U -d{U)t + v / AU-dx, (5.18) 



for any closed curve 7((U)^) that moves with the filtered fluid velocity (U)^. Thus, 
the Leray equation has an additional source of circulation arising from the difference 
between its filtered and unfiltered velocities, while the VCHE or NS-a closure model 
does not have such a term. 

This difference reappears in the vorticity dynamics for the two theories. Namely, 
upon using incompressibility V • u = we have 

<9q 

— — h u • Vq = q • Vu + i/Aq , where q = curlv, for VCHE, (5.19) 



Of 



and 



-7^- + (U)^-Voj = uj • V(U)^ + uAu) for Leray's equation 



+ 



VUjxV (U j )i , where u = curlU , (5.20) 



and we have used incompressibility of the filtered velocity, (U)^. Thus, the right 
hand side of the vorticity dynamics for Leray's equation contains an additional 
source term, compared to the curl of the VCHE or NS-a equation. 



Outlook for the remainder of the paper. We shall apply the results of Holm, 
Marsden and Ratiu Q, to recast the Lagrangian mean motion (LMM) model 
and the order 0(|£| 2 ) compressible model into the Euler-Poincare framework. Euler- 
Poincare systems are the Lagrangian version of Lie-Poisson Hamiltonian systems. 
Reformulating the LMM equations this way facilitates their Eulerian analysis, e.g., 
by providing their Kelvin-Noether circulation theorem, as well as energy and mo- 
mentum conservation as part of a general framework. We shall then use the equa- 
tions of the LMM model in making a natural adaptation and development of the 
turbulence modeling results of Chen et al. ||- Q. Namely, we shall formulate a sec- 
ond moment closure model for turbulence based on adding a certain viscosity 
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term to the ideal LMM equations. This Lagrangian mean turbulence closure model 
will then be adapted to include rotation and stratification for potential applications 
in geophysical problems. Lower-dimensional examples of Lagrangian mean theories 
will also be considered and then we shall turn our attention to developing Eulerian 
mean theories using the Euler-Poincare framework, as well. 



6 Second moment closure equations 



6.1 Euler-Poincare formulation 

We recall that the variables in the Lagrangian mean theory with averaged approx- 
imate Lagrangian (L) in equation ( |3.34j ) are the Lagrangian mean velocity u and 
the advected, or "frozen-in" quantities D and {£ k £, 1 }'- the volume element and the 
Lagrangian mean covariance of the fluctuating displacement. Such quantities sat- 
isfy a certain Lie-derivative relation Q , such as the continuity equation (|2.7| ) for the 
volume element D, 



d 

= Wt 



d 



(d 6 a) = - + £ u \ {Dd'x) = — + V-(Du) (d 6 x) , 



3D 



(6.1) 



dt u J y ' \ dt 
where £ u denotes Lie derivative with respect to the Lagrangian mean fluid velocity, 



u(x, t). There is also the geometrical relation ( 3.17 ), or its equivalent commutator 
form ( |3.25 ), for the slow time evolution of the fluctuation components for k = 
1,2,3. Namely, 



d 

- « 



.J_ 

da A 



d 
dx k 



dt + ^ * u ^)dx" 



(6.2) 



This equation for £ fc implies a geometrical relation of the same type for all the statis- 
tical moments . . . £ m ). In particular, the second moment (the Lagrangian 
mean covariance of the fluctuations) satisfies, cf. equation fl3.41| ), 



= — 



A?B\ 



(Ft 











da A da B 









dx l 



(6.3) 
(6.4) 



In vector notation this equation for the Lagrangian mean covariance dynamics is, 
cf. equation Q3.42 ), 



(6.5) 



Thus, the Taylor hypothesis ( P.17| ) provides an approximate equation for the evolu- 
tion of the fluctuation statistical moments; in particular, for their Lagrangian mean 
covariance. The Eulerian components of this equation can be rewritten as 



dt^ ? 1 



>-a t-l 



k /-a\ 



(6.6) 
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in which the right hand side is expressed as a differential operator acting on u J (x, t). 
This operator will reappear in the Lie-Poisson Hamiltonian formulation of the ideal 
second moment equations in Section [?]. Note that the isotropic, homogeneous ini- 
tial condition (£ k £ l ) = b~ kl is not invariant under the dynamics of equation ( |6,6[) 
for nontrivial velocity shear. Thus, when shear is present, the Lagrangian mean 
covariance of the fluctuations will not remain isotropic and homogeneous under the 
LMM dynamics, even if it were initially so. 

Using results of Holm, Marsden and Ratiu [jij, we compute the Euler-Poincare 
equation for the Lagrangian (L)(u, D, depending on the Lagrangian mean 

velocity u, and advected quantities D and This Euler-Poincare equation is 

given by the following extension of equation (|2.1l|), 




Thus, the additional advected quantities in general contribute their own reac- 

tive forces, appearing in the second line of the equation of motion (|6.7| ). We com- 
pute the following variational derivatives of the averaged approximate Lagrangian 
(L) in equation ( ]3.34| ) 

1 S(L) 



D 5u 

5D 
5p 



u--[d k D(ee)9i)u=v, 
-p+^i 2 + k¥>(u fc -u z ) 



l-D. 



-Pi 



tot ! 



(6.8) 



Consequently, the Lagrangian mean Euler-Poincare equation d6.7| ) for this averaged 
Lagrangian (after setting D = 1) takes the form, 



dvj 
dt 



+ U J Vi 7 j + Vju\ 



P+^|u| 2 + ^(^)(u fe -u z 



u, fc • u,) + (u,i • u,) (ee) + (u fc • u,i) (ee) 



where vi = un — (pk {^C )di) u i an d u\ = 



(6.9) 



(6.10) 



Again one sees in the entire second line of equation ( |6.S| ) the reactive forces arising 
from the variations of the Lagrangian with respect to the covariance 



Contrasting the Euler-Poincare equation ( |6.9D with the CH equation. 

In the present case, the advected second moments (£ fc £ z ) for the full system satisfy 
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their dynamical equation ( |6.5| ). If isotropic, homogeneous statistics were prescribed 
instead, so that (£, k £ l ) = a 2 5 kl with constant a 2 , then both equation (|6.5D and 



the entire second line in the motion equation (6.9) would be absent, (since the 
corresponding variations in would not be taken in this case) and the motion 

equation ( |6.S| ) of the Lagrangian mean model then would return to the n-dimensional 
generalization of the CH equation, i.e., the set ( |5,1| ) - ( |5.4| ) introduced in [[ij, Q. 
Note, however, that the n-dimensional CH equation set is not an invariant subsystem 
of the Euler-Poincare system Q6.9| ), with definition ( |6.10[ ) and advection law ( |6.5[ ), 
because (as mentioned earlier) the initial condition = 5 kl is not invariant 

under the dynamics of equation ( |6.5[ ) for nontrivial velocity shear. This means that 
the Lagrangian mean model is a departure from the VCHE model, rather than being 



an extension of it. We shall return to this matter in Section 12, when we discuss 
Eulerian mean fluid models. 

6.2 Momentum conservation — stress tensor formulation 

Noether's theorem guarantees there is a conserved momentum for the Euler-Poincare 



equations (|6.9|) , since the averaged approximate Lagrangian (L) in equation (|3.34|) 
has no explicit spatial dependence. Moreover, the integrand C in this Lagrangian is 
a polynomial in the Lagrangian mean velocity u, its gradient uj, and the advected 
quantities D and (£ k £, 1 ). That is, 

(L) = j d 3 x c(u,u k , D, , (6.11) 

with C a polynomial function of its arguments. In this case, we may express the 
Lagrangian mean Euler-Poincare equations (|6.7| ) in the momentum conservation 
form, 

d JH± = - JL T i (6 12) 

with momentum density components mi, i = 1,2,3 defined by 

_ 5(L) _dC d ( dc\ 
Su l du l dx k \ du l k J 

and stress tensor T- given by 




Equation ( |6.12| ) then implies conservation of the domain-integrated momentum, 
f md 3 x, provided the normal component of the stress tensor T- vanishes on the 
boundary. 
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In our particular case, expression (|6.14| ) for the stress tensor T- simplifies re- 
markably, to become 

2f = miu j + pSf , where m,i\ D=l = = - (d k {£, k £ l )di)ui . (6.15) 

Consequently, the equivalent Euler-Poincare motion equation (|6.9| ) simplifies to 

-i = -^ u3+p5l) - (fU6) 

The ideal Lagrangian mean motion (LMM) model. Using incompressibility 
{u? ■ = 0) reduces the expression ( |6.16| ) for momentum conservation to the La- 
grangian mean motion (LMM) equation for an ideal fluid, cf. equation ( |3.58 ), 



dv 








(6.17) 


~dt 


+ u - V v = - 


Vp, 


V • u = 0, 


where 


v - 1 W 
D 8u 


D=l 


= (i-d k {ee)9i) u. 


(6.18) 



The boundary conditions are recalled from equations ( |4.14 ), ( 4.15[ ) and ( |5.5| ) as 



v-n = , u = , and • h = 0, on a fixed boundary. (6.19) 

In ( |6.17| ) the continuity equation, V-u = 0, allows the Lagrangian mean pressure p 
to be determined from an elliptic equation, by virtue of the commutation relation 
( 3.49| ). The auxiliary equation Q6,5| ) for the slow time dynamics of the advected 



Lagrangian mean covariance may be rewritten as, 




This equation updates the Lagrangian mean covariance (£, k £, 1 ) and, thus, closes the 
LMM system. 

Contrasting the LMM model with the Euler equations. Remarkably, the 
reaction forces due to the fluctuations in the Euler-Poincare motion equation fl6.9|) 
exactly cancel the contributions of the fluctuations in the pressure and line-element 
stretching term (vjU^). This cancellation leaves only one main effect: the fluid 
parcels in the LMM model are transported by velocity u instead of v. Had we made 
no approximation at all in our Lagrangian L{uj) before averaging, we would have 
gotten exact cancellation of all of the fluctuational effects and returned entirely to 
Euler 's equations. 

So, instead of returning us entirely to the Euler equations, the Taylor series ap- 



proximation we made in Section ^3 before averaging the Lagrangian L(uj) produces 



one essential difference between the LMM model in equations ( 6. 17j ) - (6.20) and 
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the original Euler equations (2.14). Namely, the transport velocity u is smoothed 
relative to v, by inversion of the dynamical Helmholtz operator, 



u = (1 — A) v where 1 



(6.21) 



whose statistical 'metric' changes and adapts dynamically according to equa- 

tion ( |6T20D as the fluid moves. Smoothing the transport velocity to improve the 
mathematical properties of the incompressible fluid equations has been studied pre- 
viously, going back at least to the work of Leray [|^] in the 1930s, as we discussed 
earlier. However, the adaptive smoothing introduced by inverting the dynamical 
Helmholtz operator ( |6.21 ) is new here, as far as we know. Euler 's equations are 
recovered as an invariant subsystem of the LMM equation set for identically 



zero. 



We recall the remarkable commutation relation (3.49) 



d X 
dt^ D 



0, where = D~ 1 df t D(£ k £ l )di , 



(6.22) 



which may also be verified directly from the continuity equation (|2.7| ) and the La- 
grangian mean covariance equation (^). Consequently, since Ad\d = i = A, the 
ideal LMM motion equation ( 6.17] ) may be rewritten alternatively as 



du „ 
^ + U - VU 



'1 - A) _1 Vp, V 



u 



0. 



(6.23) 



Therefore, the effect of the advected fluctuations in this alternative representation of 
the LMM equation is simply to smooth the pressure gradient in an adaptive fashion 
depending on the velocity shear, through the dynamics of the covariance 

As guaranteed in advance by the mathematical theory developed in in pi], the 
Euler-Poincare equation (6.17) with definition ( |6.18| ) for the LMM model is equiv- 
alent to the Euler-Lagrange equation ( 3.44 ). Because of the commutation property 
(|6.22j ), there is also an equivalent alternative form of the LMM motion equation, 
given in ( |6.23| ). This alternative, but equivalent, form of the motion equation for 
LMM provides alternative interpretations of the effects of the fluctuation covariance 
on the Lagrangian mean motion. 



6.3 Momentum conservation for the order 0(|£| 2 ) model 
We recall the averaged approximate Lagrangian ( 3.34| ) in Eulerian form, 



(L) = |d 3 x{|[|u| 2 + (eY)(u fe -u 



+ [P 



1 d 2 p 



T - P D\. (6.24) 



2 dx k da 



Note that p is the Lagrangian mean pressure and according to equation ( 4,13j ) the 
quantity in parentheses, 



1 d 2 p if-kth 

2 dx k dx l ^ € ' 



(6.25) 
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is the Eulerian mean pressure. The momentum density components mi, i 
1,2,3 for this Lagrangian are 



mi = -^-r = D{1 - /\D)Ui , 



5u l 



where satisfies (l3.49|) and 



1 



D = i--(ee),ki 



The stress tensor T- is obtained by using equation ( 6.14 ) as 



miU 



2 dx l dx l 



(6.26) 



(6.27) 



(6.28) 



The off-diagonal components — \pn{^^ 1 ) may be regarded as a fluctuation stress 



tensor. Using the continuity equation (2/7) for D transforms the momentum con- 



servation law (6.12) into the equivalent motion equation for the LMM model with 
second order compressibility, cf. equation (|3.50|), 



(- + u-v)(l-Au>i 
with V • u 



+ 



1 dP E , 1 d ( d 2 p 

1 dD 
D~dT 



D dxi V dx l dx l 

1 d iez l ), k i + o(\^< 



2dt 



(6.29) 
(6.30) 



The boundary conditions for this model are the same as in ( |6.19| ), namely, 

v-n = , u = , and • h = 0, on a fixed boundary. (6.31) 



This completes the transformation of the Lagrangian form ( |3.44 ) of this motion 
equation to its equivalent Eulerian form. As guaranteed by the Euler-Poincare 
approach, the equivalence of equations Q3.44 ) and ( 6.29 ) may also be verified by a 
direct calculation. 



6.4 Kelvin circulation theorem for the Lagrangian mean model 



Being Euler-Poincare, the incompressible LMM equation ( |6.17| ) has a corresponding 
Kelvin- Noether circulation theorem. Namely, this equation implies 



d 
dt 



dx 



7 (u) 



7 (u) 



<9v 
~dt 



+ u • V v + VjVui 



dx 



du , 



(6.32) 



7 (u) 



for any closed curve 7(11) that moves with the Lagrangian mean fluid velocity u. This 
expression for the Kelvin-Noether property of the Lagrangian mean motion equation 
in 3D is reminiscent of corresponding expressions in wave, mean-flow interaction 
theory []l4|]. See also the Leray equation result (5.18). The main point is that 
the presence of the fluctuation covariance creates circulation of the total 

specific momentum v = u — (dk (£, k ^ l )di)u. 
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Alternatively, from equation Q6.23] ) we may write another Kelvin-Noether circu- 
lation theorem, namely 

d 



u-dxL = -f ((1- A) -1 Vp) -dx, (6.33) 

7 (u) / 7 (u) V ' 



which represents the circulation dynamics of u, rather than v. Thus, the fluctuation 
covariance creates circulation of both u and v. 

6.5 Vortex stretching equation for the Lagrangian mean model 

In three dimensions, we may use a vector identity to re-express the LMM equation 
(|6.17| ) in its equivalent "curl" form, as 

-v-ux|vxvj|Vp|u j VDj=0, V-u = 0. (6.34) 

The curl of this equation in turn yields a transport and creation equation for the 
Lagrangian mean vorticity, q = curlv, 



+ u • Vq = q • Vu + 

Ot 



V Vj x Vm j 



where q = curlv , (6.35) 



and we have used incompressibility of u. Thus, u is the transport velocity for 
the generalized vorticity q and the expected vortex stretching term q • Vu is 
accompanied by an additional vortex creation term, V vj x V u 3 . Of course, this 
additional term is also responsible for the creation of circulation of v in the Kelvin- 
Noether circulation theorem ( 6.32| ). In particular, Stokes' theorem and equation 



( 6.32| ) imply 

dS, (6.36) 



d 
dit 



curl v • dS = / / dvj Adu J 
5(u) J JS(u) 



II 


V Vj x V 


1 JS(u) 





where the curve 7(u) is the boundary of the surface S(u). One may compare this 
result with the Leray case, in equation ( 5.20| ). 



Alternatively, from equation ( |6.23| ) we may write another form of the vorticity 
dynamics, for curlu rather than curlv, namely 

^ + u-Vcj = u • Vu - V x ((1 - A^Vp) where u = curlu, (6.37) 

and we have used incompressibility of the Lagrangian mean velocity u. 

6.6 Energetics of the Lagrangian mean model 

The sum of the inner products of v with equation ( |6.23D and u with equation ( 6.34 ) 



yields conservation of energy, 

E = \f d3x {\ u I' + <^V fc • u J =^Jd 3 xu-v, (6.38) 
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after integrating by parts, using incompressibility and applying the boundary con- 



ditions ( 6.19 ). Naturally, this energy is also conserved in n dimensions. Thus, the 
covariance of the fluctuations couples to the gradients of the Lagrangian mean ve- 
locity. So it costs the system energy either to increase these gradients, or to increase 
the covariance. Moreover, there is a direct feedback between Vu and the dynamics 
of given in equation ( 6.20[). We shall see in Section ^ that Legendre transform- 
ing the Lagrangian (L) in l3~3l) gives the following Hamiltonian (still expressed 
in terms of the velocity u, instead of the momentum density m = 5{L)/5u = Dv), 



H 



d n x 



M 



D 

— I lu 

2 



+ {Z k ?)u, k -uA +p(D-l) 



(6.39) 



Remark on the geodesic property of the LMM model. When evaluated on 
the constraint manifold D = 1, the Lagrangian in ( |3.34[ ) and the Hamiltonian in 
( 6.39| ) for the Lagrangian mean motion equation ( 6.17| ) coincide in n dimensions. 
This is expected for a stationary principle giving rise to geodesic motion. The 
interpretation of the LMM model as describing geodesic motion on the volume 
preserving diffeomorphism group with respect to the H\ metric given by (U|) will 
be discussed elsewhere [Hi. 



Conservation properties of the Lagrangian mean covariance. 
equation fl6.5|) for (££), 



d 
~dl 



(^) = (a)-Vu + Vu i .(a) 



We recall 



(6.40) 



Since tr(Vu) = 0, the Lagrangian mean covariance must always have an in- 
stantaneously growing direction, along at least one principal axis of the velocity 
shear tensor, Vu. So, might be systematically growing with time. However, 
there are limits to this growth. In particular, equation ( 3.38] ) shows that this growth 
in the incompressible case must preserve the value of the determinant det along 
flow lines. Thus, stretching and rotation may occur, but for an incompressible flow 
the volume of the ellipsoid composed of the principle axes of the symmetric covari- 
ance tensor must be preserved on fluid parcels. In the Eulerian representation, 
equation ( |3.38|) implies 



d_ 

dt 



D 2 det 1=0, and 



dD 
~dt 



for V- u = . 



(6.41) 



This, in turn, implies conservation of the following quantity, for any function 



C*= [ d n x D$(D 2 det , V$, 

Jm 



(6.42) 



in which we may set D = 1 for incompressible flow. 

Conservation of the energy in equation ( 6,3§| ) also has some indications for con- 
trolling the Lagrangian mean covariance. Of course, this growth cannot continue 



indefinitely in the same direction and still satisfy conservation of energy E in (6.38) 
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in regions of nontrivial shear. So the conservation of energy must also eventu- 
ally saturate this potential growth in Lagrangian mean covariance. Thus, from the 
viewpoint of energetics, while must always have a direction in which it is grow- 
ing, this growth must also occur to preserve the determinant det and to limit 
the velocity shear Vu in accordance with the available energy in equation (|6.38| ). 
Therefore, the direction of growth (along the instantaneous principle axes of Vu 
corresponding to its positive eigenvalues) will keep changing, because of the ener- 
getic coupling and direct feedback between Vu and . A dissipative modification 
of the dynamics that allows relaxation to the homogeneous isotropic conditions 
of the VCHE model was introduced in equation ( 4.36| ). This modification may be 
used, in the event that the growth of should require additional control, e.g., in 
numerical simulations. 



7 Hamiltonian structure of the Lagrangian mean model 



Being Euler-Poincare, the LMM system consisting of the motion equation ( |6.9[) , 
the continuity equation ( |2.7[ ) and the Lagrangian mean covariance equation ( |6.6D 
must also be a Lie-Poisson Hamiltonian system. This may be verified using stan- 
dard methods, see, e.g., The corresponding Lie-Poisson bracket is dual to the 
semidirect product Lie algebra 



= u © (A e s) . 



(7.1) 



consisting of vector fields rj € u acting from the right on the direct sum of functions 
/ € A and metrics g E S. The semidirect product Lie algebra bracket has the 
expression 



[(m, h,9i),(mJ2,g2)} = ([m,mi hm - hm, aim -92m) 



(7.2) 



where we denote the induced action of u on (A © S) from the right by concatenation, 
as in /i?72- Dual coordinates are: m = 5(L)/5u = Dv dual to u; D dual to A ; 
and dual to S. The Legendre transformation of the averaged approximate 

Lagrangian (L) in ( p. 34 ) produces the Hamiltonian, 



H 



d 6 x 



^m.[D-d k D(ee)di 



-i 



m + p(D — 1) 



(7.3) 



whose proper definition requires defining the inverse of the generalized Helmholtz 
operator, (D — D(£ k £ l )di) and using the boundary conditions fl6.1S| ). The varia- 
tional derivatives of this Hamiltonian are given by 



SH 



d 6 x 



u • dm 



D 



P-2|U| 



-<£ fc £ Z }(u )fe -u,)W 



-(ufc-u,)^) + (D-l)Sp 



(7.4) 
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where evenness of the generalized Helmholtz operator under integration by parts is 
used in obtaining the first term. Finally, the Lie-Poisson bracket defined on the dual 
of the algebra in equation ( |7.1[) is given explicitly by 

{F,H}(m,D,(ee)) (7.5) 
,s f SF , n ^,5H 5F , nn .6H 5F, nn .5H 



5F 



a/-d\ rc 



5m 

5F 



5H 



In the second line of this Lie-Poisson bracket we see the contributions of the fluctu- 
ation covariance to the motion equation, cf. ( |6.7[), and in the third line we see the 
operator in the covariance dynamics itself, cf. fl6.6Q . 

The dynamical system consisting of the motion equation Q6.9D , the continuity 
equation fl2,7|) and the Lagrangian mean covariance equation ( |6.6[) now emerges as 
a Lie-Poisson Hamiltonian system in the form 



^ = {^,H}, with M G(m, £>,<£*£'» 



(7.6) 



with Hamiltonian H given in equation ( [7.3[ ) and implying conservation of the fol- 
lowing energy, cf. equation ( |3.60| ) and equation ( 6,3§| ) with m| 



E 



d 6 x 



v • ( 1 - d k (t k £ l )di 







[ d 3 x 


|u| 2 




" 2 J 







\D=1 



(7.7) 



There is a Casimir for the Lie-Poisson bracket (|7.5| ); namely, the conserved quantity 
C$ in equation ( |6.42 ). This quantity satisfies the Casimir relation, 



{C$,11} = 0, V$ and VH. 



(7.8) 



Thus, of course, C$ in equation (|6.42|) is also a constant of motion for the Lagrangian 
mean model Hamiltonian in equation (|7.3p. 



8 The 2D ideal LMM equations have no velocity Casimirs 

We recall the ideal Lagrangian mean equation of motion (|6.17j ) 

Vp, V-u = 0, (8.1) 
whore v= u- (d k (^ l )di)u= (l-A)u, (8.2) 



dv 

dt + U ' VV 



The curl form ( 6.34j ) of the motion equation is 

d_ 

dt 



■^7 v — u x (Vxvj+Vp + ^Vn^O. 



(8.3) 
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The curl of this equation, in turn, gives the vorticity dynamics, 

, where q = curlv. (8.4) 



<9q _ _ 

-f + u • Vq = q • Vu + 

at 



V Vj x Vm j 



In two dimensions these equations simplify because the vortex stretching term is 
absent, but they remain extremely nonlinear (and essentially different from the 2D 
Euler equations) because of the vorticity creation term. In two dimensions, we may 
define the stream function ?/>(x, t) satisfying u = z x V^, so that, 

q = 2 ■ q = £ • curl ((1 - A)z x W) = d x (l - A)ip x + d y (l - A)ip y = Oip . (8.5) 

This equation defines the operator O. Note that partial spatial derivatives do not 
commute with the dynamical Laplacian, A = {£, k £, l )di, because of the spatial 
dependence in (£, k £, 1 )- 

In terms of the stream function ip and the operator O, the 2D vorticity dynamics 
1) may be written as 



0_ 

Ot 



Oip + J{t/;,Oijj} = J{vj,v?} = -J{Aip x ,tp x } - J{Aip y ,ipy} , 

= -±A x \V^\ 2 y + ^A y \V^\ 2 x , (8.6) 



where J{f,g} = f x 9y — fy9x is the 2D Jacobi operation, subscripts denote partial 
derivatives and, e.g., the operator A x = {£, k £, l ) x di operates to its right. The 
highly nonlinear right hand side of this equation prevents conservation of the domain 
integrated powers of the 2D Lagrangian mean vorticity, q = Otp. For example, 
the enstrophy J q 2 dxdy is not conserved by this flow, and the Lagrangian mean 
vorticity j qdxdy is only conserved for certain boundary conditions (those for which 
§ bdy v • du = 0). Thus, the 2D ideal LMM equations apparently have no Casimirs 
that depend upon the velocity. Of course, they do have the Casimir C$ in equation 
( 6.42j ) associated with the determinant of the Lagrangian mean covariance. 



9 Relation to second moment turbulence closure mod- 
els 

By adapting the work of Chen et al. [Q]- ||, we form a second moment turbulence 
closure model that combines the one point closure model they studied with the re- 
sults of the previous sections and yields the following system of equations, 



— v + (u-V)v = i/Av-Vp + F, V-u = 0, (9.1) 



where v is a constant kinematic viscosity and in this case v is given by 

v = u - (d k (£ k t l )di)u = (1 - A)u . (9.2) 
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Note that the Lagrangian mean fluctuation covariance appears in the dissipation 
operator A. In the absence of the forcing F, this viscous LMM turbulence model 
dissipates the energy E in equation ( |6.38 ) according to 



dE 
~dt 



j d 3 x tr(Vu T • • Vu) + Au • Au . (9.3) 



This strictly negative energy dissipation law is the reason for adding viscosity with 
A, instead of using the ordinary Laplacian operator. 

These equations are closed by the dissipative dynamical equation ( |4.36| ) for the 
advected Lagrangian mean covariance, rewritten as, 

^<^> = u k 3 (ee) + ^ - - «eY> - ^ + a— a(^) . (9.4) 

at ' J ,J T t 

The boundary conditions for this dissipative LMM model for a fixed boundary 
are v = 0, u = 0, for the velocities and either • h = when dissipation is 
absent in the dynamics of or equation ( p.4| ) without the diffusion term when 

such dissipation is present. 

As mentioned before, isotropy and homogeneity of the Lagrangian mean co- 
variance tensor are not preserved under the dynamics of equation (|9.4| ) for 
nontrivial velocity shear. Therefore, the VCHE or NS-a closure model of Chen et 
is not an invariant subsystem of the dissipative LMM model in equations 
).4[) . However, the phenomenological addition of relaxation and dissipation to 
the covariance dynamics does make the dissipative LMM model relax to the VCHE 
or NS-a closure model, whose validity in steady state has been verified by compar- 
ison with experimental data, ||— ||]. These dissipative LMM equations comprise 
a 3D dynamically self-consistent second-moment turbulence closure model, whose 
solution may be sought as a dynamical systems problem, e.g., as an initial value 
problem for decaying turbulence, or as motion of a forced dissipative system. 

The symmetric quantity n$ = d(££) jdt in equation fl9,4|) is called the Taylor 




diffusivity [18] and is sometimes used in turbulence modeling as a semi-empirical 
model for passive scalar diffusion (T^j , [|l9| , |5| . This application will be considered 
further in Section [H], when we include the effects of rotation and stratification, as 
well. For now, the main effect of the fluctuations in equations (|9.1| ) - ( |9.2[) from the 
viewpoint of the Navier-Stokes equations for an incompressible fluid is to smooth 
the transport velocity u relative to the circulation velocity v, by inversion of the 
Helmholtz operator (1 — A). This smoothing tends to suppress triad interactions at 
wave numbers greater in magnitude than about _1 / 2 . Therefore, the cascade 

of the spectral kinetic energy density E(k) = ±(u7v)(k) as a function of wave 
number will be suppressed for |k| > |(££)| -1 / 2 . Hence, viscous dissipation should 
take over at length scales smaller than the local length scale given by | | 1//2 (x, t) 
and suppress the effects of nonlinearity at these smaller scales, just as if the LMM 
model in (^TTJ) - Q were an adaptive LES scheme. The LES aspects of the 
dissipative LMM model will be pursued elsewhere. 
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10 Geophysical applications — adding rotation and strat- 
ification to the Lagrangian mean model 

10.1 Formulation of Lagrangian mean Euler-Boussinesq (LMEB) 
equations 

At leading order in |£| 2 , introduction of rotation and stratification alters the aver- 
aged approximate Lagrangian (L) in equation ( p. 34 ) to the following form, 



(L) 



d A x{ — 



D r 



|u| 2 + (^)(u fc .u/)J +DR(x)-u-gbDz+ p 
The corresponding Euler-Poincare equations are, cf. [0, ||, 



1 - D 



(10.1) 



d 
dt 



v — u x curlR(x) + Vp + gbz = , with V- u = , 



(10.2) 



where — = (— + u -v), v = (l-A)u, A = V • • V , 
|=0, and ^{«) = <«)-Vu + Vu T .(«), (10.3) 



with boundary conditions, 



v-n = u = , and • n = on a fixed boundary. 



(10.4) 



These Lagrangian mean Euler-Boussinesq (LMEB) equations describe the 
Lagrangian mean effects of fluctuations on the ideal motion of a stratified incom- 
pressible fluid in a rotating reference frame, given in the Eulerian description. (Note 



that the equation in ( [L0.3D does not refer to internal waves. See [14| for a dis- 
cussion of the Lagrangian mean effects of internal waves from a viewpoint similar 
to the one taken here.) The motion equation for this system may be rewritten as 



d 

— v — u x curl (v + R) + u^'Vvj + Vp + gbVz 



0. 



(10.5) 



From this expression, we immediately find the Kelvin- Noether circulation the- 
orem showing how the covariance of the fluctuations interacts with the buoyancy 
to generate mean total circulation. Namely, 



d_ 

dt 



(v + R) • dx 



7(u) 



5(u) 



V Vj x V u 3 - gVfexV z 



dS, 



(10.6) 



where the closed curve 7(11) moves with the Lagrangian mean fluid velocity u and 
is the boundary of the surface S(u). The terms [V«j XV u 3 ] and [g V6 xVz] 
are both vorticity creation terms in this Lagrangian mean Kelvin-Noether cir- 
culation theorem. Therefore, the fluctuations act together with buoyancy to create 
circulation of the sum v + R around closed fluid loops moving with velocity u. In 
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particular, the vorticity creation term [V Vj x V u 3 ] may cause mixing across isopy- 
cnal (constant b) surfaces, even when the buoyancy is stably stratified. According 
to the interpretation given in Section ||, this is the effect of the Stokes mean drift 
velocity. 

The fluctuations and the buoyancy also act together in the conserved energy 
for the LMEB system (JT02|) - ( pip . Namely, cf. equation ( jTllD , 



Elmeb 



d 6 x 



u| 2 +tr(Vu T -(£0 • Vu) +gbz 



(10.7) 



This expression shows the trade offs in energetics among translational kinetic energy, 
Lagrangian mean velocity shear combined with the fluctuations, and gravitational 
potential energy that may occur in the dynamics of the LMEB model. These en- 
ergetic tradeoffs should have implications for the stability properties of the LMEB 
model's equilibrium solutions. We note that the LMEB model has all the same 
equilibrium solutions of the original Euler-Boussinesq equations (namely, those on 
the invariant manifold = 0) in addition to others associated with critical points 
of the sum of the energy Elmeb and the conserved Casimirs for the EAB theory, 

Clmeb = J d 3 x D$(b, Ab, D 2 det (££)) for an arbitrary function, <£> . (10.8) 

These quantities Poisson commute with every Hamiltonian H: that is, under a 
modification of the Lie-Poisson bracket in equation (ffl]) to include the buoyancy 
variable 6, as follows, 



{F,H} LMEB = {F,H}^- J d' 



x 



5F . , S SH S F / \ 5H 



5nii 5b 5b 

we have {C LM eb, H}lmeb = 0, for all H[m,D, (££},&]. 



5m a 



(10.9) 



10.2 Passive scalar diffusion in the LMEB model. 



Following Bennett [19], we introduce the symmetric Taylor diffusivity tensor [18] as 



2dt^ W 



2dt 



(10.10) 



According to the semi-empirical theory [T 
scalar equation with this diffusivity is 

db 







- = V-(K S .V b). 



, the corresponding passive 



(10.11) 



As indicated here, we shall adopt this equation for the dissipative dynamics of 
the buoyancy, even though b is not strictly passive. The present theory provides 
a dynamical equation ( |10.3| ) for the evolution of the Taylor diffusivity tensor, Kg- 
Thus, passive scalars in this theory advect and diffuse according to equation ( |10.11| ) 
with tensor diffusivity given by equation (|6.40| ) as 



2k s = j t m = <«> • Vu + Vu T • («> 



(10.12) 
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Hence, the passive scalar equation we find is, 

1/ V-((^)-Vu + Vu T -{^))-V)/). (10.13) 



alt 2 

where the Lagrangian mean covariance and the velocity shear tensor Vu are 
determined self-consistently from the LMEB dynamics, suitably modified to include 
dissipation. We see that the tensor diffusivity k,s = " ^ u + ^ uT, (£0) 

is determined dynamically as the symmetrized product of the Lagrangian mean 
covariance and the mean velocity shear of the flow. As we have discussed, equation 
( 10.12 ) implies that the covariance at any time will be increasing along the 
instantaneous positive eigendirections of the velocity shear tensor Vu. Conservation 
of energy ( 10. 7|) should control this tendency, in general, because the cost in energy 
for the covariance to grow is highest precisely where the growth rate is highest. 

The covariance dynamics following from the Taylor hypothesis for the advection 
of turbulent structures by the (Lagrangian) mean flow implies the determinant of 
wu l be conserved along flow lines in the incompressible case. Therefore, the 
covariance will remain finite and nonzero. However, if any additional dissipation is 
needed to moderate the growth of due to shear forcing, equation (0) may be 
used as an alternative expression for the Taylor diffusivity 

Finally, following Chen et al. Q- ||], we introduce viscosity into the motion 
equation for the Lagrangian mean Euler-Boussinesq theory as A-diffusion of mo- 
mentum v, cf. equation (|9.1D, with boundary conditions as discussed after equation 
(|9.4| ). This final step - the semi-empirical introduction of dissipation - completes 
our formulation of the LMEB equation set ( |1 . 1| ) - ( |1.6| ) for the Lagrangian mean 
motion of a rotating stratified incompressible fluid. In what follows, we shall con- 
sider some lower dimensional subcases of this Lagrangian mean motion in Section 
11, then turn our attention to the corresponding Eulerian mean theory in Sections 
12 and O. 



Remark on asymmetry of the diffusivity and other modeling steps. Be- 
cause of the Earth's rotation, asymmetric diffusivity may be expected in the ocean. 
This is an additional modeling step, which will be pursued elsewhere. For dis- 
cussions of this aspect, see, e.g., the review by Davis and references therein. 
Other modeling steps would address the creation of fluctuation covariance, as well 
as its dissipation, or treat the covariance between fluctuations at two different times. 
We do not pursue these other sub-gridscale modeling directions here. Instead, we 
proceed from the Taylor hypothesis ( |3.17| ) that the Lagrangian mean fluctuation 
covariance is advected by the Lagrangian mean fluid velocity. 



11 Lagrangian mean motion in fewer dimensions 

11.1 2D Lagrangian mean rotating shallow water (LMRSW) equa- 
tions 

To the averaged approximated LMEB Lagrangian (L) in equation ( |10.1| ) we ap- 
ply the standard shallow water approximations: neglect vertical gradients; neglect 
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kinetic energy of vertical motion; assume constant density; and integrate in the 
vertical direction from the bottom topography at z = —B(x,y) to the free sur- 
face at z = h(x,y,t). These approximations lead to a new Lagrangian (L)^ sw for 
Lagrangian mean rotating shallow water (LMRSW) dynamics, 



(L) L RS W = J dx dy | ^ [|u | 2 + <£Y) (u fc • u , 



+ DR(yL)-u-±gD 2 + gDB{x,y)\. (11.1) 

Here D = h + B(x,y) is the total depth of the water. The corresponding Euler- 
Poincare motion equation for LMRSW dynamics in momentum conservation form 
is, cf. equations ( |6.12j ) - ( |6.14[ ) , 

= - -^j( m ^ + \9D 2 5l) +gDB >i (x,y)+Du j B j ti , (11.2) 
where one defines momentum density components m^, for i = 1, 2 by 

m, L = ^^f 1 = D{ Ui + Ri) - {D(^ k ^) Uiyl ) k = D(l - A D )u t + DRi . (11.3) 
The depth D satisfies the continuity equation, 

^ + V-( J Du) = 0, (11.4) 

and (^^) satisfies the Lagrangian mean covariance dynamics inherited from the 
Taylor hypothesis ( |3.17j ), 

|<«> = <«>-Vu + Vu T .<€0 • (11.5) 

In terms of the circulation velocity v = (1 — Ad)u defined as before, but now in 
2D, we find the motion equation for v, 

— + u •v)v-uxcurlR(x)+ 5 v( J D-S(a;,y)) =0. (11.6) 

This is the same as the standard motion equation for rotating ideal shallow water 
dynamics, modulo the substitution u — > v and the additional dynamics for the 
covariance This motion equation for LMRSW implies the Kelvin-Noether 

circulation theorem, cf. equation ( 6.32; ). 

d_ 

dt 



(v + R)-dx= / / dvjAdv? = / / Vti 3 xVii 3 -(iS. (11.7) 

7(u) J Js(u) 



If 


V Vj x V11 3 


1 JS(u) 





Thus, as for the LMM model, the fluctuations have the effect of generating (total) 
circulation. The curl of the LMRSW motion equation yields, with q = z-curl (v+R), 

(l + U • V )(l) =^ l£ -( V ^i x V « =D- 1 J{v i ,v?). (11.8) 
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Thus, perhaps unexpectedly, the potential vorticity q/D is not conserved on fluid 
parcels by the LMRSW model, but instead has a local creation term, D~ 1 j(vj,u^. 
This creation of potential vorticity is in fact the convective mechanism by which the 
smoothing of velocity gradients occurs in the LMRSW model. The usual equations of 
rotating shallow water theory are recovered on the invariant subsystem for = 0, 
on which v = u, and the right hand side of the potential vorticity equation Qll.8 ) 
vanishes. 



11.2 One dimensional Lagrangian mean shallow water (LMSW) 
equations 



Restricting the Lagrangian (L)^ sw in equation ( 11.1 ) to one dimensional motion 
(without rotation) results in 



(L) 



sw 



dx 



D I 2 2 
— [U +WU X 



-gD 2 + gDB(x) 



(11.9) 



Here we denote by w in one dimension. We also denote partial derivatives by 
subscripts and D = h + B(x) is the total depth of the water. The correspond- 
ing Euler-Poincare equations in momentum conservation form for LMSW in one 
dimension are, cf. equations (|6.12|) - ( |6. 14 ) , 



m t 

with m defined by m 
where D and w satisfy D t 
and wt 



mu H — q D' 



+ gDB x , (11.10) 



6(L) L 



SW 



5u 
-{Du) x , 

-u w x + 2w u x 



Du — (Dwu x ) x . 



(11.11) 
(11.12) 
(11.13) 



The corresponding motion equation for the momentum velocity v = m/D is given 

by, 



vt + uv x + g[D — B(x)) x = , where 



u 



D l (Dwu x ) x . 



(11.14) 



Thus, the Lagrangian mean equation in the case of one dimensional shallow water 
again makes just one simple change v — > u in the advection term, and introduces an 
additional equation for the dynamics of the Lagrangian mean covariance (££) = w. 

An explicit calculation using equations ( 11.1 2| ) , ( |11,13| ) and the definition of the 
momentum velocity v in ( |11.14 ) regains the remarkable relation ( 0.49 ), expressed in 
one dimension as, 



v t + uv x = (1 - D 1 d x Dwd x )(u t + uu :L 



IH.15) 



Hence, the one dimensional Lagrangian mean shallow water equation system may 
be rewritten as 



u t + uu x + g (1 - D- l d x Dwd x )~ l {D - B) 
D t + uD x + Du x = 0, 
wt + u w x — 2w u x = . 



0, 



(11.16) 
(11.17) 
(11.18) 
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This system for w 7^ is no longer hyperbolic; rather, it is conservative and disper- 
sive. However, when w = 0, we again recover the usual shallow water dynamics as 
an invariant subsystem. 

One may also obtain the equations for Lagrangian mean polytropic gas 
dynamics with pressure-density relation p = po(D / Z?o) 7 , by replacing the motion 
equation in the Lagrangian mean shallow water system ( |11.16| ) - ( 11.1 8| ) with 



Ut+UU, 



+ M(i 
+ D? [ 



D-^drDwdrT^LV-^D, = 0. 



(11.19) 







Thus, equations ( |11.17|) - ( 11.19|) provide a system of Lagrangian mean poly- 
tropic gas equations, which is no longer hyperbolic for w 7^ 0, but which recovers 
the usual polytropic gas dynamics as an invariant subsystem when uu = 0. This 
system conserves the energy, 



E 



dx 



2 2 
u + wu„ I + 



Po 



jyi 



7 



(11.20) 



The one dimensional Lagrangian mean models for shallow water and polytropic 
gas dynamics given in this Section should provide ample opportunities for testing 
the dynamical effects of fluctuations on one dimensional Lagrangian mean fluid 
motion, without requiring the more intensive numerics and analysis needed in higher 
dimensions. (The three dimensional version of the Lagrangian mean polytropic gas 
equations is readily obtained at this point.) 



11.3 The Lagrangian mean Riemann (LMR) equation 

The process we apply here of fast-slow decomposition, followed by Taylor approxi- 
mation and averaging at fixed Lagrangian fluid label in Hamilton's principle for ideal 
fluids produces the following situation. There are more equations in the averaged 
approximate model than in the original model and they include the original model 
as an invariant subsystem in which the fluctuation covariance vanishes. 

The simplest situation for which this occurs is probably the Riemann equation, 

u t + 3uu x = 0, (11-21) 

which describes one dimensional shock formation and is the Euler-Poincare equation 
for the Lagrangian, 

L R = j dx X -u 2 . (11.22) 

The corresponding averaged approximate Lagrangian for this problem is (again writ- 
ing = w for the covariance in one dimension) 

(L) R = ! dx ^(u 2 + wu 2 x ), (11.23) 
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whose Euler-Poincare equation is 



1 2 1 2 
v t +[-u + -wu x + uv 



with 



6(L) 



R 



5u 



u - {wu x ) x , 



and the one dimensional Lagrangian mean covariance w satisfies 

wt = — u w x + 2w u x . 



(11.24) 



(11.25) 



(11.26) 



The Lagrangian mean Riemann system (11.24) - (11.26) can be rewritten equiva- 
lently in (nonlocal) characteristic form as 



u t + uu x 
and wt 



l-d x w d x ^j d x (u 2 - -w u x ! . 
uw x + 2w u x ■ 



(11.27) 
Jll.28) 



We compare the Lagrangian mean Riemann system (11.24) - (11.26) with the 
Camassa-Holm (CH) equation, a one dimensional completely integrable nonlinear 
dispersive shallow water model [21], [22|, 

(u - u xx ) t + {^u 2 - ^u 2 x + u(u- Uvx)} = , (11.29) 

which may also be written in (nonlocal) characteristic form closely similar to equa- 
tion ( |11.271) as 



u t + u u x 



l-dl 



d x (u 2 + ^u 2 x 



(11.30) 



The Camassa-Holm equation ( 11.29| ) is the Euler-Poincare equation for the La- 
grangian 



f 1 2 

Lch = I dx -(u + u 



(11.31) 



Further discussions of the solution behavior for the Lagrangian mean Riemann sys- 
tem (11.24) - ( 11.26j ), including its traveling waves (which may be solved explicitly) 
will be given elsewhere [42|. For now, we simply observe that there is a key difference 
in sign between the LMR equation (11.27) and the CH equation ([L1.30]). 



12 Eulerian mean theory of advected fluctuations 

Here we develop an alternative Euler-Poincare theory of advected fluctuations from 
the viewpoint of Eulerian averaging. We have the same point of departure as for the 
Lagrangian mean theory, namely the Lagrangian ( 2.10| ) in the Eulerian description, 



L(u) 



yd 3 x|^|U(x,t; W )| 2 +P(x,t; W )(l- J D(x,f;u;))| . (12.1) 
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The traditional Reynolds decomposition of fluid velocity is expressed at a given 
position x in terms of the Eulerian mean fluid velocity, U as 

U(x, t;u) = U(x, t) + U'(x, t;u). (12.2) 

According to equation ( [4.6D the Eulerian velocity fluctuation U'(x, t;u) is related 
to the Eulerian displacement fluctuation — denoted as £(x, t ; u) — by 



dC - 

-± + U- VC = C • VU + U'(x, t ; uj) 



(12.3) 



As discussed in Section ^ for purely Eulerian velocity fluctuations as in equation 
( 12.2j ), this relation separates into two relations: the "Taylor- like" hypothesis of [23], 



and the relation 



§+u.v C =o; 



= C- VU + U'(x,t;w). 



(12.4) 



(12.5) 



Hence, the Reynolds velocity decomposition ( |12.2 ) separates the Lagrangian ( |12.1[ ) 
into its mean and fluctuating pieces, as 



L(u) 



d 3 x { — |U(x, t) + U'l 2 + P(x, t)[l- D(x, t 



(12.6) 



No modification is needed in the pressure constraint in this Lagrangian, because the 
Eulerian mean preserves the condition that the velocity be diver genceless; hence, 
V-U = 0. It remains only to take the Eulerian mean of this Lagrangian, in which 
we assume (C) E = 0. The Eulerian mean averaging process at fixed position 
x is denoted ( • ) E with, e.g., 

i r T 

U(x, t) = (U(x, t ; uj)) e = lim - / U(x, t ; u) du . (12.7) 



By equation ( |12.5 ), the Eulerian mean kinetic energy due to the velocity fluctuation 
satisfies 



IU 



/|2 \£ 



(c fc cr u k ■ u, z . 



(12. 



Thus, we find, the following Eulerian mean Lagrangian, (cf. equation ( 3.54| ) for the 
corresponding Lagrangian mean form) 



d*x{ — 



D r 



|U(x, t)| 2 + ( C\ l ) E ^,k ■ U,i + P(x, t) { 1 - £>(x, t) 



(12.9) 



The advection relation ( p.2.4 ) implies a similar advection of each component of the 
symmetric Eulerian mean covariance tensor {( k ( l } E - Th us, we have 



^ + u-v)<cW = o. 



(12.10) 
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This relation the continuity equation for the volume element D, 

dD 

— + V-DXJ = 0, (12.11) 

complete the auxiliary equations needed for deriving the equation of motion for the 
Eulerian mean velocity tj from the averaged Lagrangian (L) E in ( [12.9 ) by using the 
Euler-Poincare theory. 

The results of Holm, Marsden and Ratiu allow one to compute the Euler- 
Poincare equation for the Lagrangian (L) (tj, D, (C k ( l ) ) m ( |12.9| ) depending on 
the Eulerian mean velocity tj, and advected quantities D and (( k ( l ) E as, cf. equa- 
tions (pUP and flBTTp, 



d -, d \ 1 5(L) E 1 S(L) E -,• 

TZ + U's-t) 7^477r + T^^TtV^ (12-12) 



dt dxi J D 5U l D SW 

_d_ 5^f 1 8(L) E _d_ k l E 
dx i 5D D 5{( k ( l ) E dx i ^ ^ ; 

kA\E 



Thus, the Eulerian mean covariance ^ k Q) h ' saLisfying Lhe scalar adveclion rela- 



tion ( 12. 10|) c ontributes considerably simpler reactive forces than those arising in 



equation (|6.7| ) from the Lagrangian mean covariance satisfying the tensor advection 



equation (pj). We compute the following variational derivatives of the averaged 



approximate Lagrangian (L) in equation ( 12.9 ) 



1 5(L) E _ T 1 (a „,, k , hE: 



u-- [d k D(c K Crd l u = v 



D 5U D 

5 -^f = -^ + ^|u| 2 + ^(CW(u,-uO--^ 

6(L) E 



5P 

5(L) E D 



l-D, 

S(CH l ) E §(U,fc-U,0 ■ (12.13) 

Of course, these variational derivatives are in the same form as in equation set ( |6.8[ ). 
However, the variables here are Eulerian mean quantities and they will enter a dif- 
ferent Euler-Poincare equation, namely, ( 12.12| ), instead of ( |6.7j) . This difference 



arises because the Eulerian mean covariance (C C / advects as an array of scalars 
under the Eulerian mean evolution, rather than as the components of a symmetric 
tensor. Consequently, the Euler-Poincare equation ( |12.12| ) for this averaged La- 
grangian takes the form, 



^ + WV hJ + Vjlji = - - \ (U * ■ tjji) (C k C l )° , (12.14) 

1 



01 ' ' ' dx 

~D 



where V l = U i - hf,U l with = j- (d k D (( k ( l ) E di) and 1^ = 0. (12.15) 



The boundary conditions we shall choose for this motion equation are 

V-n = , U = , and h • (££} E = 0, on a fixed boundary. (12.16) 
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Then, provided the Helmholtz operator 1 — for D = 1 may be inverted, the 
Eulerian mean pressure P may be obtained by solving an elliptic equation. 



Contrasting the Euler-Poincare equation ( 12.14 ) with the CH equation. 

When the Eulerian mean covariance is isotropic and homogeneous, so that (C k C l ) E = 
a 2 5 kl , for a constant length scale a, then this equation reduces to the n-dimensional 
Camassa-Holm equation introduced in mi, S. Thus, the n-dimensional CH equation 



set is an invariant subsystem of the Euler-Poincare system ( 12.14 ), with defini- 
tion (12.15) and advection law ( 12.10| ), because the initial condition (C k C l ) E = cx 2 5 kl 
is invariant under the dynamics of equation ( |12.1C| ). 

Physical interpretation of V as the Lagrangian mean velocity for EMM. 

(12.17) 



The Stokes mean drift velocity is defined by [12], 

(U) s e((.VU') £ . 
Hence, equation ( |12.5| ) implies 

(U) 5 = - (C • VC ■ V) S U = - A E \J + o(\C\ 2 ) 

where 

A E = (d k (C k ( l ) E di) = Af 



;i2.i8) 



(12.19) 



and we again argue that V • C = °(IC| 2 )- Thus, we find that V satisfies, to order 
o(IC| 2 ), 



V EE U 



•-u = U + (U) 



(U) 



(12.20) 



Therefore, to this order, V is the Lagrangian mean velocity for the EMM 
theory. Thus, the duality between the Lagrangian mean velocity and the Eulerian 
mean velocity is reciprocated. In Eulerian mean theories the dual momentum is the 
Lagrangian mean velocity, and vice versa. 



Contrasting the Euler-Poincare equation (12.14) with the CL equation. 



The EMM motion equation ( 12.14 ) may be rewritten equivalently as 



V V+ (V-U)xcurlV + VP 

v * ' 

Stokes vortex force 



i<CV) B V(U fe 



Ui) =0. 



(12.21) 



Thus, the Stokes mean drift velocity V — tj = A E JJ contributes an additional 
vortex force in this motion equation for the Lagrangian mean fluid velocity V. 
This equation is similar in form to the CL equation ( |4.26| ), However, the meaning 
of the velocities are reversed in the two cases: the CL equation is for the Eulerian 
mean velocity U; and the EMM equation is for the Lagrangian mean velocity V. 
The Stokes mean drift velocity is also prescribed for the CL equation, rather than 
being determined dynamically. 
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12.1 Kelvin circulation theorem for the Eulerian mean model 



Being Euler-Poincare, the Eulerian mean model (EMM) equation ( |12.14 ) has a cor- 



responding Kelvin- Noether circulation theorem. Namely, this equation implies, 
cf. equation (|6.32|), 



'' 1 V-dx = -i / / _ VfUfe-U,) x V(C k C l ) E -dS, (12.22) 



dt / 7 (tj) 2 j js(u) 

for any closed curve 7(U) that moves with the Eulerian mean fluid velocity U 
and surface S(Xj) with boundary 7(U). Thus, in this Kelvin-Noether circulation 
theorem the presence of the Eulerian mean fluctuation covariance {C k C l ) E creates 
circulation of the Lagrangian mean velocity V = U - A E XJ. 

12.2 Vortex stretching equation for the Eulerian mean model 

In three dimensions, the EMM equation ( 12. 14| ) may be expressed in its equivalent 
"curl" form, as 

(9„ T - T /_ „,\ fT _ \rt _ 1 /tt . tt \T7lAkA\E 



-V-Ux (Vx VJ +V (P^ + U-V) = - -(u fc -Ui)V(CTr, v-u = o. 

(12.23) 

The curl of this equation in turn yields an equation of transport and creation for 
the Lagrangian mean vorticity, Q = curl V, 

^+U- VQ = Q- VU- iv(U fe .Uj) X V(( k C l ) E , where Q^curlV, 

(12.24) 

and we have used incompressibility of U. Thus, U is the transport velocity for 
the generalized vorticity Q and the expected vortex stretching term Q • VU is 
accompanied by an additional vortex creation term. Of course, this additional 
term is also responsible for the creation of circulation of V in the Kelvin-Noether 



circulation theorem ( 12.22 ) and it vanishes when the Eulerian mean covariance is 
homogeneous in space, thereby recovering the corresponding result for the three 
dimensional CH equation Q, @|. 

12.3 Energetics of the Eulerian mean model 

Noether's theorem guarantees conservation of energy for the Euler-Poincare 
equations (12.14), since the Eulerian mean Lagrangian {L) E in equation ( |12.9| ) has 



no explicit dependence on time. This constant energy is given by 

E = \f d 3 x(|U| 2 + (C fc C Z ) £ U fc .U z ) =\ y>xU-V. (12.25) 

Thus, the total kinetic energy is the integrated product of the Eulerian mean and 
Lagrangian mean velocities. In this kinetic energy, the Eulerian mean covariance 
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of the fluctuations couples to the gradients of the Eulerian mean velocity. So there 
is a cost in kinetic energy for the system either to increase these gradients, or to 
increase the Eulerian mean covariance. 

As one might expect from the analysis in Section fj], Legendre transforming the 
Lagrangian (L) E in (|12.9| ) gives the following Hamiltonian (still expressed in terms 
of the velocity tj, instead of the momentum density M = S(L) E /5XJ = DV), 



H 1 



d n x 



M 



+ ((Y) £ U r U i )+P(fl-l) 



(12.26) 



Remark on the geodesic property of the EMM model. When evaluated on 
the constraint manifold D = 1, the Lagrangian in (12.9) and the Hamiltonian in 
(12.26) for the Eulerian mean fluid motion equation ( 12.14]) coincide in n dimen- 



sions. This is expected for a stationary principle giving rise to geodesic motion. 
The interpretation of the EMM model as describing geodesic motion on the volume 
preserving diffeomorphism group with respect to the H\ metric given by ( |12.25| ) is 
discussed in |23| for the CH case, in which the Eulerian mean covariance is isotropic 
and homogeneous. The corresponding discussion for the EMM equations is equiv- 
alent to that given in [23], since the metrics are equivalent, provided the initial 
conditions for (( k C l ) E are bounded away from zero. 

12.4 Momentum conservation — stress tensor formulation 

Noether's theorem also guarantees conservation of momentum for the Euler-Poincare 



equations ( 12.14 ), since the Eulerian mean Lagrangian {L) E in equation ( |12.9| ) has 
no explicit spatial dependence. As before, the integrand C in this Lagrangian is a 
polynomial in the Lagrangian mean velocity U, its gradient U fc, and the advected 
quantities D and {( k ( l ) E ■ That is, 

(L) E = J d 3 xC(tJ,lJ <k , A<CW), (12-27) 
with C a polynomial function of its arguments. In this case, we may express the 



Eulerian mean Euler-Poincare equations ( 12.14| ) in the momentum conservation 
form, 

dM % d 



dt dxi 1 ' 

with momentum density components Mi, i = 1,2,3 defined by 

5{L) E dC d ( dC \ 



Mi 



5U l 



dW dx k \ dm 



(12.28) 



(12.29) 



and stress tensor given by, cf . equation Q6.14 ) , 
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Equation ( 12.28| ) then implies conservation of the total momentum, f'M.d^x, pro- 
vided the normal component of the stress tensor TP vanishes on the boundary. 
In our particular case, expression (12.30) for the stress tensor TP becomes 



TP = MiW - DtJ, h • U i (C k ( j ) E + Pdf , where Mi = D{U % - K E Ui) . (12.31) 



Consequently, the equivalent Euler-Poincare motion equation ( |12.14| ), or ( |12.23| ) is 
also expressible as 



dt 







Vi W + P5{ - Ujfc • U ti (CV) E 



0riV , , } . with Vi = Mi\ D=x . (12.32) 

The boundary conditions are given in equation ( 12.16j ). 



12.5 A second moment turbulence closure model for EMM 



When dissipation and forcing are added to the EMM motion equation (12.14) by 
using the phenomenological viscosity vK E "V and forcing F, one finds a second 
moment Eulerian mean turbulence model given by 



d 



1 



kA\E 



_ + tJ.V)V + ^W + VP t %+ -(U fc -U,i)V<C fc CT (12.33) 



i/A^V + F, where V-U = 0, 



with viscous boundary conditions V = 0, U = at a fixed boundary. Note that the 
Eulerian mean fluctuation covariance (C k C^) E appears in the dissipation operator 
A E . In the absence of the forcing F satisfying appropriate regularity conditions, 
this viscous EMM turbulence model dissipates the energy E in equation fljggg ) 
according to 



dE I ^3 

— = — v I d 6 x 

dt 



tr(VU J • (CC) -VU) + A fc U • A^U 



(12.34) 



This negative definite energy dissipation law justifies adding viscosity with A E , 
instead of using the ordinary Laplacian operator. 



13 Geophysical applications of the Eulerian mean model 
13.1 Eulerian mean Euler-Boussinesq equations 

Introducting rotation and stratification alters the averaged approximate Lagrangian 
(L) E in equation ( 12. 9| ) to the following expression, 



(L) 



E 

EB 



d A X 



D 
~2L 



|u|" + <c fe C 



kA\E 



U k • u , 



+ DR(x) -U -gbDz+ P 



1 - D 



(13.1) 
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To account for buoyancy b in the Euler-Poincare equation ( 12.12j ) , the right hand side 
should have the additional summand, D~ 1 b t i5(L) BB /5b. Then, the Euler-Poincare 
equation resulting from (L) BB in ( |1 3.1 ) is 



<9V 

— + U-VV + VjVW 
at 



gbz - VP; 



U x curlR(x 
1 



tot 



kA\E 



where V = {l-K E )\], with A £ = V • ( CO -V 



(13.2) 
(13.3) 



and the Eulerian mean flow is incompressible, so V-U = 0. The auxiliary equations 
that complete this set are 



0, (! + U-V)(CO' 



3D 

0, ¥ + V-(DU) = 0, (13.4) 



These Eulerian mean Euler-Boussinesq (EMEB) equations describe the 
Eulerian mean effects of fluctuations on the ideal motion of a stratified incompress- 
ible fluid in a rotating reference frame. The motion equation for this system may 
be rewritten in curl form, as 



3 V - U x curl(V + R) + V(P t E ot + U • V) + gbVz + l - (U fc ■ U,) V^C^ 



dt 



0. 



(13.5) 



The Kelvin- Noether circulation theorem for the EMEB model is given by, 

1 



d_ 

1ft 



(V + R) • dx 



7(U) 



S(U) 



V(U fc .Uj) xV(CVr+ gVbxV 



dS. 
13.6) 



Here, the closed curve 7(U) moves with the Eulerian mean fluid velocity U and is 
the boundary of the surface S(U). There are two vorticity creation terms in the 
Kelvin-Noether circulation theorem ( 13.6 ). Thus, spatial variation in the Eulerian 
mean covariance and nonvertical buoyancy gradients can both create circulation 
of the sum V + R around closed fluid loops moving with velocity U. Naturally, 
the same vorticity creation terms appear in the dynamics of the Lagrangian mean 
vorticity Q = curl V, which may be obtained by taking the curl of equation (|13.5| ). 
The conserved energy for the EMEB system ( (DTp - (fOOD 

is given by a 

formula similar to equation ( 10.7| ) for the LMEB model, which is 



Eemeb 



d 3 x 



i(|u| 2 + tr(vu T .(cc) E -vu; 



+ gbz 



(13.7) 



Adding dissipation semiempirically as before yields equations ( 1.18Q - ( |1.19| ) in the 
Introduction. 
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13.2 2D Eulerian mean rotating shallow water (EMRSW) 



We apply the standard shallow water approximations to the averaged approximated 
EMEB Lagrangian (L) E in equation ( 13.1| ) to find a new Lagrangian (L)^ sw for 
the dynamics of Eulerian mean rotating shallow water (EMRSW), 



(L)rsw 



dxdy ^ — ^ 



|u|^ + (cV) E (u fe -u 



1 



+ DR(x)-\J-^gD 2 + gDB(x,y)j. (13.8) 

Here D = h + B{x, y) is the total depth of the water, which satisfies the continuity 
equation, 



dD 
~dt 



+ V-(DU) = 0, 



(13.9) 



and the Eulerian mean covariance (££) E satisfies the dynamics inherited from the 
Taylor-like hypothesis ( 12.4j ), 



(JUo.v)<CC> E = 0. 



;i3.io) 



The Euler-Poincare motion equation ( 12.12| ) for EMRSW dynamics generated by 
the Lagrangian {L)^ sw is expressed as 



|+o. V IV 



U x curlR(x) + gV[D- B{x,y) 



VjVU 3 



-(U fc .U,0V(C fe Cr , (13.11) 



where the Eulerian mean circulation velocity V = (1 — A^)U is defined as in 
equation ( 12.15| ), but now in 2D. The left hand side of this equation is the same 
as the standard motion equation for rotating ideal shallow water dynamics, modulo 
the substitution tj — > V and the additional dynamics for the covariance (CO^"- 
The right hand side has a "line element stretching term" and a "covariance gradient 
term." 

This motion equation for EMRSW implies the Kelvin- Noether circulation 
theorem, cf. equation ( |12.22j ), 



^ <b _ (V + R) • dx 



dt 



7 (U) 



5(U) 



V(U fe -U z ) xV(C fc C 



kA\E 



dS. 



[13.12) 



with fluid loop 7(U) and surface S(\J), as before. Thus, as for the EMM model, the 
covariance gradient term has the effect of generating (total) circulation. The curl of 
the EMRSW motion equation yields, with Q = £-curl(V + R), 



\dt 



"+0.V)(§ 



1 



^s.v(u, fc .u,,) xv(c k c l ) E 



;i3.i3) 
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Thus, the potential vorticity Q/D is not conserved on fluid parcels by the EMRSW 
model, but instead has a local creation term proportional to V(C fc C') S - However, 
conservation of circulation on fluid loops and conservation of potential vorticity Q /D 
on fluid parcels is recovered for the invariant subsystem ( CC) E = o?5 kl with constant 
a 2 . For this invariant subsystem, we have V = (1 — a 2 A)U (with the ordinary 
Laplacian operator) and the covariance gradient ^((^ k ^ l ) E vanishes. Hence, the 
right hand sides of equations ( |13.12 ) and ( |13,13| ) also vanish in this case; so that 
the Kelvin circulation integrals are constant and the potential vorticity is conserved 
on fluid parcels in this case. 



13.3 ID Eulerian mean shallow water model 

Restricting the Lagrangian (L)^ sw in equation ( |13.8j ) to one dimensional motion 
(without rotation) results in 



(L) 



E 

SW 



dx 



D 



u 2 + wu: 



-gD 2 +gDB(x) 



(13.14) 



Here W denotes ( CC) E i n ID, subscripts denote partial derivatives and D = h+B(x) 
is the depth of the water. The corresponding Euler-Poincare equation is, from 
equation ( 12.12| ), 



v t + uv x + vu x + -w x u 2 + 



g(D — B{x)) - -(U 2 + WU : 



(13.15) 



This may also be written in momentum conservation form for EMSW in ID as, cf. 



equations ( p^8|) - ( p^0|) 

M t = 



(Mil -DWU 2 + -gD 2 ) + (,DB.,.. 



with M defined by M 

where D and W satisfy Dt 
and Wt 



SW 



6U_ 
- (DU) X , 

-uw x . 



Dtj - (DWU X ) X 



(13.16) 

(13.17) 

(13.18) 
(13.19) 



The equivalent equations (13.15) and ( |13.16| ) simplify into the following motion 
equation for the Lagrangian mean velocity V = M/D, 



V t + UV x - D~ l (DWU 2 



where V = U 



-g(D-B(x)) x 
D- l (DWU x ) x . 



(13.20) 



Thus, the Eulerian mean equation for one dimensional shallow water makes the 
change V — » U in the advection term, and introduces another term on the left hand 
side. Amusingly, this other term is half the commutator of the advective derivative 
d/dt and the Helmholtz operator 1 — acting on U. That is, 



^ , ( w ®-[(l +[? 3'( 1 - fl t D ^ 



d 







u . 



(13.21) 
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Therefore, we have two equivalent forms of the ID EMSW equations 
V t + UV x -D-\DWU 2 x ) x = -g{D-B(x)) 



(1 - Ag)(0t + UU X ) + D-\DWU 2 



x Jx 
x)x 



X 



(13.22) 
(13.23) 

where V = (1 - A£) U = U - D~\DWU X ) X . (13.24) 
Hence, the one dimensional Eulerian mean shallow water system may be 



-g(D-B(x)) 

-l 



Ut 


+ uu x 


+ (1 - Ag)- 1 [g (D - B) x + D-\DWU 2 X ) X ] = , 


(13.25) 


A 


+ UD 2 


+ DU x = 0, W t + UW x = 0, 


(13.26) 


(i 


-Ag) 


= (l-D-^DWd*). 


(13.27) 



This system for W ^ is no longer hyperbolic; rather, it is conservative and dis- 
persive. When W is constant, we obtain the shallow water alpha-model as an 
invariant subsystem. 

The equations for Eulerian mean polytropic gas dynamics with pressure- 
density relation p = po(D / Dq)" 1 , are obtained by replacing the motion equation in 
the Eulerian mean shallow water system ( [13.25] ) - ( 13.27] ) with 



U t + UU x + 1 



P^D^ 2 D X + D- 1 (DWU 2 ) X 



0. 



(13.28) 



Thus, equations ( [13.26] ) - (13.28) provide a system of Eulerian mean polytropic 
gas equations, which is no longer hyperbolic for W ^ 0, and which recovers the 
polytropic gas alpha model as an invariant subsystem when W is constant. This 
system conserves the energy, 



E 



dx 



7 



(13.29) 



13.4 The Eulerian mean Riemann (EMR) equation 



For another illustration of the Eulerian mean methodology in one dimension, we 
return to the Riemann equation, 



U t + 3UU X = 0. 
This is the Euler-Poincare equation for the Lagrangian, 



Lr 



dx -U 2 



(13.30) 



(13.31) 



The corresponding Eulerian mean approximate Lagrangian for this problem is (again 
writing ( CC) E = W for the Eulerian mean covariance in one dimension) 



(L)j 



dx l -(JJ 2 + WU 2 x ), 



(13.32) 
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whose Euler-Poincare equation is 



V t + UV X + 2VU X + ~W X Ui = . 



with 



V 



5_{L)_ E 
5U 



R 



u-(wu x ) x , 



(13.33) 



(13.34) 



where the one dimensional Eulerian mean covariance W satisfies the scalar advection 
law 



W t 



uw x 



(13.35) 



The Eulerian mean Riemann system (13.33) - ( |13.35| ) can be rewritten in momen- 
tum conservation form as 



Vl + \2 U 



W -n - 

-jU 2 x + UV 



0. 



and in nonlocal characteristic form as 

Ut + UUx = 
and Wt = 



l-d x Wd x 
UW X . 



(13.36) 



(13.37) 
(13.38) 



We compare the Eulerian mean Riemann system ( 13.33| ) - ( 13.35 ) with the com- 
pletely integrable model of Camassa-Holm |2^| for nonlinearly dispersive shal- 
low water waves. The equation of motion for the CH model can be written in 
momentum conservation form as, cf. equation ( 13.36 ), 



(u - u xx ) t + (\u 2 - \u 2 x + U{U- U xx ) 



0. 



(13.39) 



Moreover, the CH equation may also be written in nearly the same (nonlocal) char- 
acteristic form as equation ( |13.37| ) for the EMR problem, 



u t + uu x 



l-tt 



\ (u 2 + \m 



(13.40) 



The CH equation ( |13.3S| ) is the Euler-Poincare equation for the Lagrangian 



Lch 



dx-(U 2 + Ul) 



(13.41) 



Thus, the Eulerian mean Riemann system ( 13.33 ) - ( |13.35 ) reduces to the CH equa- 
tion ( |13.3S| ) for W = 1 (or for any other nonzero constant, which can be absorbed 
into the spatial length scale). Hence, the Eulerian mean Riemann system is the 
natural extension of the CH model to allow for time dependence of the length scale 
associated with W, the Eulerian mean covariance of the rapid fluctuations. More 
analysis of this equation is given in 42]. 
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14 Conclusions 

We have considered two classes of models that describe the mean motion of a fluid 
in the presence of rapid or random fluctuations. These are: 

• The Lagrangian mean fluid models, in which the fluctuation is modeled 
as a displacement of the Lagrangian fluid parcel trajectory and averages are 
taken holding its Lagrangian label fixed; and 

• The Eulerian mean fluid models, in which the fluctuation is modeled as 
occurring at a fixed position, based on the traditional Reynolds decomposition 
of the fluid velocity, and averages are taken at fixed spatial position. 

At linear order in a Taylor expansion in the magnitude of the fluctuations we obtain 
the relation that allows us to treat the two classes of models on the same 

mathematical footing. We apply asymptotics and averaging methods to Hamilton's 
principle for an ideal fluid and use the Euler-Poincare theory of Holm, Marsden and 
Ratiu jij], |2] to determine the equations of motion that result in each class of model. 

From the Euler-Poincare viewpoint, the two classes of models differ from each 
other primarily in the way they treat Taylor's hypothesis, that the fluctuations 
should be "frozen" into the mean flow. We take Taylor's hypothesis to mean that 
the fluctuations should "transform" in a certain way under the action of the mean 
flow. In the Eulerian mean models, the displacement fluctuation transforms under 
the action of the Eulerian mean flow as a collection of scalars, by parallel transport. 



(See [23 1 for a mathematical description of this action using parallel transport in 
a composition of diffeomorphisms in a manifold setting.) In the Lagrangian mean 
models, the displacement fluctuation transforms under the action of the Lagrangian 



mean flow as a vector field, by Lie transport as in equation ( 3.17| ). 



The two different models yield different equations via the Euler-Poincare theory. 
However, these equations are "dual" to each other in the physical interpretations of 
their solutions. Namely, the momentum evolving in the Lagrangian mean models 
is interpreted in Section |I| as the Eulerian mean velocity. And, vice versa, the 
momentum evolving in the Eulerian mean models is interpreted in Section ^ as 
the Lagrangian mean velocity. The energy in both models is the total mean kinetic 
energy. The formal source of this duality between the two classes of models turns 
out to be their shared form of conserved total mean kinetic energy. This kinetic 
energy in both cases may be interpreted as the domain-integrated product of the 
Eulerian mean velocity times the Lagrangian mean velocity. These two velocities 
are related in both classes of models by a dynamical Helmholtz operator whose 
metric is the covariance of the fluctuations. Evenness of this Helmholtz operator 
allows the kinetic energy to be written as the H 1 norm of whichever velocity is being 
studied. Then, the variational derivative of the kinetic energy with respect to one 
of these velocities summons the other one and thereby produces the duality, in their 
Euler-Poincare equations. 

The effect of the averaging in either case is to make the solution velocity smoother 
than the momentum, or circulation velocity that it transports, via the inversion of 
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the Helmholtz operator that relates the two velocities. The covariance that provides 
the metric appearing in this Helmholtz operator introduces a length scale which 
evolves with the mean flow according to the appropriate Taylor hypothesis for the 
fluctuations. The magnitude of this length scale determines the smoothness of the 
solution velocity for both classes of models. The mechanism for the smoothing 
is nonlinear dispersion in these ideal fluid theories (before viscosity is added). 
The nonlinear dispersion contains the length scale associated with the covariance 
and acts to suppress the magnitude of the fluidic triad interaction at smaller length 
scales. Normally, the triad interaction in fluid dynamics drives the cascade of energy 
forward to smaller length scales. However, in both of the classes of mean fluid models 
we study here this process is suppressed by nonlinear dispersion. In summary, 

The nonlinear dispersion in these mean fluid theories acts to make the 
transport velocity smoother than the circulation or momentum velocity, 
by the inversion of a Helmholtz operator whose length scale corresponds 
to the fluctuation covariance. 

These mean theories each provide either an extension, or a development of the 
viscous Camassa-Holm equation (VCHE, or NS-a) that has recently been introduced 
as a one-point turbulence closure model, Q- ||. The Eulerian mean models are 
indeed natural extensions of the VCHE, or NS-a models to second order closures for 
turbulence. The Eulerian mean models reduce to the VCHE, or NS-a model when 
the Eulerian mean fluctuation covariance is spatially homogeneous. The Lagrangian 
mean models are another departure which also provides a second moment closure for 
turbulence that is related to the VCHE, or NS-a model, but does not contain it as 
an invariant subsystem unless the covariance vanishes entirely. Which approach will 
eventually lead to an appropriate model for climate and other long time geophysical 
applications remains to be seen. 

We also formulated several examples of these Eulerian mean and Lagrangian 
mean fluid models in fewer dimensions, in the hopes that these simpler examples 
will provide additional physical insight into the dynamical behavior of these two 
classes of models. 
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